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INT  RODL’CT  ION 

A  periodic  structure  is  one  which  consists  cf  a  number  cf  identical  substructures,  called 
periodic  units,  coupled  together  to  from  the  entire  structure.  There  are  in  general  two  types  cf 
periodic  structures  -  simple  periodic  structures  and  compound  periodic  structures.  A  simple  periodic 
structure  is  one  which  consists  cf  basic  periodic  units  which  cannot  be  divided  further  into  identical 
subunits.  Figure  1  shows  an  example  cf  such  a  structure.  Figure  1  shows  a  uniform  beam 
supported  at  spacing  t.  Oscillators  are  attached  at  raids  pan  between  the  supports  as  shown.  The 
oscillators  are  modeled  as  consisting  cf  mass  (md),  elastic  stiffness  (spring  constant  kd)  and 
damping  (dashpot  constant  cd)  A  compound  periodic  structure  consists  cf  periodic  units  which  in 
themselves  are  composed  cf  identical  periodic  subunits.  Thus  a  periodic  unit  in  a  compound  periodic 
structure  is  itself  another  periodic  structure.  Figure  2  shows  an  example  cf  a  compound  periodic 
structure.  Figure  2  shows  a  uniform  beam  simply  supported  at  spacing  t  where  oscillators  are 
equally  spaced  within  the  spans  as  shown. 

Periodic  structures  have  been  analyzed  using  the  concepts  cf  a  'propagation  constant" 
[1-3]  and  also  by  the  '’transfer  matrix  method"  [4],  The  concept  cf  propagation  constants  as  applied 
to  wave  propagation  in  periodic  structures  is  reviewed  in  Appendix  A.  However,  the  derivation  of 
propagation  constants  is  very  cumbersome,  requiring  the  solution  cf  forth  •  cxder  partial  differential 
equations  in  most  cases.  Cn  the  ether  hand,  due  to  recent  advances  in  transfer  matrices  [5,5],  the 
transfer  matrix  method  appears  to  be  less  cumbersome.  This  approach  permits  a  simple  treatment 
d'  periodic  units  with  complicated  configurations  and  furthermore,  a  matrix  formulation  1  .cost 
suitable  for  periodic  structures  cf  finite  total  length  since  the  imposition  cf  boundary  conditions  at 
the  ends  cf  the  structure  is  straight  forward  [7]. 

In  this  report,  an  attempt  is  made  to  utilize  the  transfer  matrix  method  for  analyzing 
wn'e  propagation  in  periodic  structures.  Three  examples  are  given  to  illustrate  the  general 
ippraich.  The  cise  of  longitudinal  vibration  and  wave  propagation  in  a  rod  is  considered  first.  Thas 
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is  fallowed  by  mare  complicated  problems  which  include  vibration  and  wave  propagation  in  both  a 
two-dimensional  3-bay  planar  lattice  structure  and  a  three-dimensional  tetrahedral  truss. 

AN.\LY  SIS 

TOE  TRANSFER  MATRIX  METHOD 

State  Vectors 

The  identification  of  a  state  vector  is  important  in  applying  the  transfer  matrix  method 
to  wave  propagation  and  vibration  of  structures.  The  state  vector  z  at  a  point  i  of  an  elastic  system 
is  a  column  vector,  the  components  ctf  which  are  the  displacements  at  the  point  and  the 
corresponding  internal  forces. 

For  example,  in  the  longitudinal  vibration  ctf  a  straight  rod,  the  state  vector  z  consists  of 
components  u  and  ,V ,  where  u  is  the  longitudinal  displacement  and  N  is  the  axiai  fcrce.  As 
another  example,  for  the  analysis  of  the  flexural  vibration  of  a  Timoshenko  beam,  the  state  vector 
consists  ctf  components  w ,  M  and  V  where  w  is  the  transverse  displacement,  if/  is  the  rotation  Ln 
radians  ctf  the  cress-section,  and  M  and  V  are  the  moment  and  the  shear  force,  respectively. 

In  the  analysis  ctf  a  structure,  if  the  time  histories  ctf  the  state  vectors  for  specified 
locations  in  the  structure  are  known,  the  vibration  characteristics,  as  well  as  the  wave  propagation 
characteristics,  can  be  determined. 


Transfer  Matrices 


A  transfer  matrix  relates  the  state  vector  at  a  point  in  a  structure  to  the  state  vector  at 
another  print  in  the  same  structure.  It  is  an  nxn  matrix,  where  n  refers  to  the  number  ctf 
components  Ln  the  corresponding  state  vector. 


Some  common  transfer  matrices  are  derived  in  .Appendices  B  through  E.  The  transfer 
matrices  .'from  left  to  rich:  and  from  rich:  to  left,  the  siariticance  cf  which  is  explained  later  Ln  this 
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section)  far  longitudinal  vibration  in  an  elastic  rod  with  distributed  mass  are  derived  in  Appendix  B. 
The  transfer  matrix  for  longitudinal  vibration  in  an  elastic  rod  with  distributed  mass  and  damping  is 
derived  in  Appendix  C.  The  transfer  matrix  for  flexural  vibration  in  an  elastic  beam  with 
distributed  mass,  including  the  effects  of  shear  deflection  and  rotary  inertia  is  derived  in  Appendix  D 
and  the  transfer  matrix  far  both  longitudinal  and  flexural  vibrations  in  a  bar  with  distributed  mass, 
including  the  effects  of  shear  deflection  and  rotary  inertia,  is  derived  in  Appendix  E. 

For  example,  the  transfer  matrix  for  longitudinal  vibration  in  a  straight  rod  relates  the 
state  vectors  at  the  two  end  points  as  (refer  to  Appendix  B) 


cos  9 


sin9 

9 


IsinQ 
EA  9 
cos  9 


(1) 


and  for  flexural  vibration  in  a  Timoshenko  beam,  the  transfer  matrix  relates  the  state  vectors  at  the 
two  end  prints  as  (refer  to  Appendix  D) 
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where  the  subscripts  R  and  L  denote  the  right  and  left  end  state  variables  faces,  respectively.  For 


-14- 


tte  definitions  of  the  variables  used  in  eqns.  (1)  and  (2),  refer  to  Appendices  B  and  D  respectively. 

Appendix  C  also  shows  the  general  approach  for  deriving  the  transfer  matrix  for 
longitudinal  vibration  in  an  elastic  rod  with  damping.  However,  since  the  introduction  c£  damping 
complicates  the  transfer  matrices  considerably,  the  problem  cf  damping  will  not  be  considered  here. 

Uses  of  Transfer  Matrices 

As  mentioned  previously,  a  transfer  matrix  is  iced  to  relate  the  state  vectors  at  two 
specified  points  in  a  structure.  Notice  that  due  to  the  sign  conventions  chosen  in  deriving  transfer 
matrices,  a  transfer  matrix  becomes  different  when  a  left  end  state  vector  cf  a  particular  element  is 
written  in  terms  cf  the  right  end  state  vector  as  opposed  to  when  the  right  end  state  vector  is  written 
in  terms  cf  the  left  end  state  vector. 

For  example,  for  longitudinal  vibration  in  a  straight  rod  (refer  to  Appendix  B),  the 
transfer  matrix  from  left  to  right  relates  the  state  vectors  at  the  two  ends  cf  the  rod  as 


cos  9 


l  sin  9 
EA  9 
cos  9 


and  the  transfer  matrix  from  right  to  left  relates  the  state  vectors  at  the  two  ends  as 


(3) 


ccs  9 
sin  9 


9 


_  l  sin  9 
EA  9 
ccs  9 


It  can  be  shewn  that  the  transfer  matrix  in  eqn.  (4)  can  also  be  obtained  by  simply  taking  the 
inverse  of  the  transfer  matrix  in  eqn.  (3)  and  vise  versa  (refer  to  eqn.  (B16)). 


One  cf  the  useful  and  interesting  characteristics  cf  the  transfer  matrix  method  is  that 
transfer  matrices  can  be  multiplied  to  form  another  transfer  matrix  which  represents  a  larger  section 
cf  a  structure.  .As  a  simple  example,  consider  the  red  cf  length  2 i  os  shewn  in  Figure  3.  The  rod  is 


assumed  to  consist  cf  two  identical  elements,  each  cf  length  t.  Using  eqn.  (3), 
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and 
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Substituting  eqn.  (5)  into  eqn.  (6)  gists 


l  sin  9 

cos  0  EA  9 
sin  9  cos  0 


cos  0 


l  sin  9 
EA  0 
cos  9 


(7) 


Since  A  and  C  are  the  left  and  right  end  points  cf  a  rod  cf  length  It,  eqn.  (7)  states 
that  the  transfer  matrix  for  longitudinal  vibration  cf  a  rod  cf  length  It  is  equal  to  the  product  of 
two  transfer  matrices  for  a  rod  cf  length  t.  To  prove  the  validity  cf  such  a  statement,  first  calculate 
the  product  of  the  two  transfer  matrices.  Some  mathematical  manipulations  give 
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ccs  0  EA  9 

cos  29  “  EA  29 

.  *  sin  9  ccs  9 
^  e  ] 

p  i  sin  9  ccs  9 
.  9 

2^  sin  29  055  29 

-0 

(8) 


Since  9  —  t  ■Vi  (refer  to  Appendix  B) , 


29  =  2fu>\/  -2-,  eqn.  (7)  can  be  obtained  by  substituting  It  for  t  in  the  transfer  matrix  in  eqn. 


-16- 


(3).  Thus,  the  transfer  matrix  for  longitudinal  vibration  cf  a  rod  of  length  21  is  indeed  equal  to  the 
product  cf  two  transfer  matrices  for  a  rod  of  length  l. 

In  fact,  the  technique  cf  transfer  matrix  multiplication  can  be  applied  to  more 
complicated  structures  since  the  intermediate  state  vectcrs  can  be  substituted  successively  to  obtain 
the  transfer  matrix  for  the  entire  structure.  It  is  this  characteristic  which  makes  the  transfer  matrix 
method  a  favorable  approach  in  analyzing  periodic  structures.  Thus,  transfer  matrices  can  be 
combined  in  such  a  way  that  intermediate  state  vectors  can  be  eliminated. 

For  a  specific  problem,  the  excitation  and  the  boundary  conditions  must  be  specified. 
When  an  excitation  is  applied  to  a  structure,  say  at  point  (or  station)  p,  the  state  vector  becomes 
discontinuous  at  p.  The  problem  is  solved  using  the  boundary  conditions.  This  is  illustrated  in  the 
next  paragraph.  For  example,  consider  the  rod  shown  in  Fig.  4.  The  rod  consists  cf  four  identical 
rod  segments  cf  length  i  and  is  loaded  by  a  sinusoidal  axial  force  of  magnitude  Na  at  point  C. 

Now  assume  that  the  rod  in  Fig.  4  is  broken  at  C .  Fig.  5  shows  the  forces  at  the  left 
end  and  the  right  end  cf  the  rod  at  C .  Since  displacements  are  continuous  in  crossing  point  C  (that 
is,  uL  at  C  =  u„  at  C), 


where  Cm  and  CL  denote  points  just  to  the  right  and  just  to  the  left  cf  C,  respectively. 

Let  T  (<)  be  the  transfer  matrix  represented  by  the  2x2  matrices  in  eqn.  (3)  when  i  in 
parentheses  signifies  that  T  is  for  the  transfer  matrix  cf  a  rod  cf  length  i.  At  CL , 


=  7- (07(0 


(10) 


Since  the  rod  is  assumed  to  be  a  continuous  member  with  no  impedance  mismatch,  multiplying 


T  (€)  n  times  is  equivalent  to  replacing  (  by  n(  in  T .  Thus,  eqn.  (10)  can  be  rewritten  as 


/“l  M 


Using  eqn.  (11),  eqn.  (9)  becomes 


However,  for  point  E , 


Canbining  eqn.  (12)  and  eqn.  (13)  gives 


fcl  -feo[reofc}4  ♦  {&_ 


u\  fO  \ 

'  v  '  =  W  v  +  T  (2€) '  v  (14) 

K‘y)t  vY  Ja  {‘yo  Jc 

at  the  boundaries,  since  the  the  displacements  are  specified  (that  is,  u  =0),  eqn.  (14)  can  be 


wntten  as 


(°L  +7-™{U 


Eqn.  (15)  can  be  solved  to  obtain  the  internal  forces  at  ^4  and  E .  In  ether  words,  the 
state  vectors  at  A  and  E  can  be  obtained  by  imposing  the  boundary  conditions. 
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Since  this  is  steady  state  vibration,  the  state  vector  at  intermediate  points  in  the  rod  can 
be  obtaired  using  either  are  c£  the  boundary  state  vectors  [8]  (that  is,  the  state  vectors  at  A  or  £). 
For  example ,  for  point  B , 


(16) 


or 


and  far  point  D , 


=  T~1(3() 


+  T(t) 


(17) 


(18) 


or 


(19) 


where 


r-‘(0 


_  ^  sin  9 
cos  9  £.4  9 

„  1  sin  9  cos  9 

— — 


is  the  transfer  matrix  in  eqn.  (3). 

In  addition,  frequency  response  functions  for  ether  specific  locations  in  a  structure  can 
be  obtained  through  the  use  c£  transfer  matrices  [8].  With  the  frequency  response  functions  known, 
random  vibration  can  be  considered.  Moreover,  impulse  response  functions  can  be  generated  from 
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tte  frequency  response  functions,  whereby  time  histories  cf  waves  in  a  structure  can  be  studied.  [10] 

The  relationships  between  transfer  matrices  and  propagation  constants  in  periodic 
structures  has  been  investigated  in  [9].  It  has  been  found  that  at  any  particular  frequency,  the 
propagation  constants  corresponding  to  the  waves  in  a  structure  are  equal  to  the  negative  natural 
logarithms  cf  the  eigenvalues  cf  the  transfer  matrix  relating  the  state  vectors  at  the  two  ends  cf  the 
basic  dement  constituting  the  periodic  structure.  Since  propagation  constants  give  information  on 
attenuations,  wave  numbers  and  phase  changes  for  wave  propagation  in  structures  (refer  to 
Appendix  A),  the  wave  propagation  characteristics  can  be  readily  obtained,  via  the  propagation 
constant  technique  once  the  transfer  matrix  far  a  periodic  unit  is  derived. 
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APPLICATION  OF  THE  TRANSFER  MATRIX  METHOD  TO  AN  ELASTIC  ROD  WITH  DISTRIBUTED 
MASS 

Figure  6  shows  the  rod  to  he  analyzed  in  this  section.  The  rod  has  modulus  c i  elasticity 
£ ,  mass  density  p  and  crass-sectional  area  A .  In  addition,  the  rod  is  assumed  to  be  made  up  of  six 
identical  rod  segments,  each  cf  length  £;  thus  there  are  no  impedance  mismatches  throughout  the 
length  of  the  rod. 

State  Vectors 

For  longitudinal  vibration  in  a  red,  the  state  vector  consists  cf  a  longitudinal 
displacement  component  u  and  an  axial  farce  component  N .  Figure  7  shows  the  sign  convention  to 
be  used  in  this  analysis.  Thus, 


(20) 


Transfer  Matrices 

The  transfer  matrices  (from  left  to  right  and  from  right  to  left)  for  longitudinal 
vibration  in  an  elastic  rod  with  distributed  mass  are  derived  in  Appendix  B.  From  eqn.  (Bll),  the 
transfer  matrix  T  which  relates  the  state  vector  at  the  right  end  to  the  state  vector  at  the  left  end  cf 
an  elastic  rod  cf  length  t  is  given  as 


T  - 


— uf 


ccs  9 
sin  9 
9 


i  sin  9 
E.\  9 

cos  9 


Cl) 


and  from  eqn.  (316),  the  transfer  matrix  7_1  which  relates  the  state  vector  at  the  left  end  to  the 
state  vector  at  the  right  end  cf  an  elastic  rod  is  given  as 
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f  sin  9 

cos  0  EA  0 
p£(j  ^  Q  cos  0 

where  the  variables  in  eqn.  (21)  and  (22)  are  defined  in  Appendix  C. 

Since  the  rod  to  be  analyzed  in  this  section  has  a  length  cf  6f ,  and  since  the  rod  is 
assumed  to  have  constant  material  and  geometric  properties  throughout  its  length  with  no 
impedance  mismatches,  it  is  convenient  to  divide  the  rod  into  six  periodic  units,  with  each  periodic 
unit  represented  by  a  rod  cf  length  l.  Thus  eqns.  (21)  and  (22)  become  the  transfer  matrices  (from 
left  to  right  and  from  right  to  left)  cf  one  periodic  unit  cf  the  system. 

Uses  of  the  Transfer  Matrix 

With  the  transfer  matrices  far  one  periodic  unit  cf  the  rod  defined,  the  state  vector  far 
specific  locations  cf  the  rod  can  be  obtained.  Far  example,  referring  to  Fig.  6,  the  state  vector  at  C 
going  from  left  to  right  along  the  rod  is  given  by 


a: 
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where  the  values  in  parentheses  indicate  the  appropriate  arguments  in  terms  cf  t  in  the  matrix  T. 
Notice  that  if  there  are  impedance  mismatches  in  the  rod,  no  such  simplification  can  be  made.  For 
example,  if  section  AB  cf  the  rod  is  made  up  erf  a  material  with  modulus  cf  elasticity  2 E  instead  cf 
£  as  for  the  rest  of  the  rod,  eqn.  (23)  becomes 


where  T  represents  the  transfer  matrix  for  section  AB  and  is  obtained  by  substituting  2£  for  £  in 
matrix  T .  Notice  that  eqn.  (24)  is  not  affected  since  section  AB  is  not  included. 

Dumping  can  also  be  included  in  the  transfer  matrix  method.  The  transfer  matrix  from 
left  to  right  far  longitudinal  vibration  in  an  elastic  rod  with  distributed  mass  and  damping  is  derived 
in  Appendix  C.  Assume  that  the  rod  in  Figure  6  now  has  material  damping  which  can  be 
characterized  by  a  viscous  damping  constant  c  and  has  no  impedance  mismatches.  The  state  vector 
at  C  can  still  be  obtained  using  eqn.  (23)  tail  with  the  transfer  matrix  T  defined  by  eqn.  (C  11)  of 
Appendix  C  (instead  of  eqn.  (22)). 

The  frequency  response  functions  for  specific  locations  in  the  rod  are  derived  in 
Appendix  F.  Appendix  F  also  contains  the  general  approach  whereby  the  problem  cf  forced 
vibration  in  a  rod  can  be  treated.  The  rod  is  assumed  to  be  excited  by  a  sinusoidal  axial  force  cf 
magnitude  .V,  at  E  as  shown  in  Figure  3.  Fran  eqn.  (F  16), 


=  - 


cos  28 
cos  60 


where  the  subscript  A  denotes  the  response  location  and  the  superscript  £  denotes  the  exdtaticn 
location  and  the  subscript  N  denotes  an  axial  fcrce  response.  Frcm  eqn.  (F22), 


-t 
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The  impulse  response  function  at  B  is  also  generated  in  Appendix  F.  Due  to  the  simple 
nature  cf  tff( ui),  the  impulse  response  function  is  obtained  in  dosed  form.  From  eqn.  (F  44),  the 
impulse  response  for  the  force  at  B  is  given  as 

**1(0  =y£  {-^K3+24*)fVJj  -^r-(5+24*KVf] 

-*f  -(7+24 k  )l VJ]-B(f  -(9+24 ■*)<  VJ] 

-tf[r  -<15+24jkK\/l']+«(f-(17+24Jfc)€V^|] 

E  E 

+6 (t  -<19^24  k  -(21  +24* 

E  E 

(36) 

For  an  impulse  excitation  at  E,  eqn.  (36)  gives  the  time"  history  of  the  response  to  be 
expected  at  B.  Two  points  are  cf  particular  interest  here.  First,  the  fraction  y  in  front  cf  the 

summation  sign  signifies  that  the  amplitudes  cf  the  responses  at  B  with  respect  to  time  is  always  one 
half  of  a  delta  function.  This  is  expected  because,  due  to  symmetry  in  the  rod,  the  axial  force 
excitation  is  split  up  into  two  equal  and  opposite  going  waves  traveling  along  the  rod.  Second, 

notice  that  the  terra  has  the  unit  cf  time.  In  fact,  it  represents  the  time  required  for  the 

axiai  wave  to  travel  a  distance  l  along  the  rod. 


In  addition,  at  any  particular  frequency,  the  propagation  constant  can  be  obtained  as  the 
negative  logarithms  cf  the  eigen  values  cf  transfer  matrix  7.  Since  7  is  a  2x2  matrix,  there  is  only 
one  pair  cf  equal  and  opposite  propagation  constants,  corresponding  to  two  opposite  and  identical 
waves  as  noted  earlier  [9]. 
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APPLI  CATION  OF  THE  TRANSFER  MATRIX  METHOD  TO  A  3- BAY  PLANAR  LATTICE  STRUCTURE 

Egure  9  shows  the  3-bay  planar  lattice  structure  to  be  analyzed  in  this  section.  For 
simplicity,  assume  that  the  cross-sectional  dimensions  in  the  bars  are  small  compared  with  the 
lengths,  and  that  the  structure  is  made  up  of  identical  horizontal  and  vertical  bar  elements 
throughout.  Furthermore,  the  Timoshenko  beam  model  is  used  fcr  the  bars  such  that  the  effects  ctf 
shear  deflection  and  rotary  inertia  are  included  in  the  analysis. 

State  Vector 


Recall  that  the  vibration  and  wave  propagation  in  a  structure  are  characterized  by  a 
state  vector  z.  In  the  case  of  the  planar  structure  shown  in  Figure  9,  for  each  bar,  t»  e  state  vector  z 
consists  of  three  displacement  component  and  three  internal  force  components.  The  three 
displacert^nt  components  are  u ,  w  and  ill,  where  u  is  the  longitudinal  displacement,  w  is  the 
transverse  displacement  and  vj»  is  the  rotation  of  the  cross-section.  The  three  force  components  are 
M ,  V  and  .V,  where  M  is  the  moment,  V  is  the  shear  farce  and  IV  is  the  axial  force.  Figure  10 
shows  the  sign  convention  for  the  forces  and  the  displacements.  Thus 


,  =  d 

‘  V 


where  d  is  the  displacement  vector  and  p  is  the  force  vector  such  that 


(37) 


and 


P  = 


In  analyzing  the  structure  shown  in  Figure  9,  state  vectors  corresponding  to  both  main 
members  I  and  n  are  needed  to  describe  the  vibration  characteristics  of  the  structure.  The  reason 
for  this  is  due  to  the  choice  ctf  the  transfer  matrices  and  will  become  apparent  in  the  next  section. 
Thus,  for  a  particular  section  in  the  structure,  the  state  vectcr  Z  cf  interest  is 
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Z  = 


d„ 

Pu 

Pi 


(38) 


where  the  subscripts  denote  member  numbers.  Notice  that  Z  is  used  to  differentiate  the  state  vector 
at  a  section  in  a  structure  from  that  c£  z,  which  represents  the  state  vector  at  a  point  in  the 
structure. 


Transfer  Matrices 

Two  transfer  matrices  are  involved  here.  The  first  transfer  matrix  Xu  involves  the 
transfer  of  state  vectors  in  two  tars,  each  of  length  l  in  main  members  I  and  n.  For  example, 
members  12  and  12'  together  are  represented  by  such  a  matrix.  The  second  transfer  matrix  X2 
involves  the  transfer  cf  state  vectors  across  the  junctions.  The  members  which  join  main  members  I 
and  II  constitute  such  a  matrix.  For  example,  member  11  constitutes  such  a  matrix.  Figure  11 
shows  the  3-bay  planar  structure  which  has  been  sectioned  into  its  constituent  pars  responsible  for 
transfer  matrices  X,  and  X2.  The  subscripts  R  and  L  are  used  to  denote  prints  which  are  just  to 
the  right  and  points  which  are  just  to  the  left,  respectively,  cf  junctions  which  join  main  members  I 
and  II. 


Transfer  matrices  X;  and  X2  are  derived  in  Appendix  G.  From  eqn.  (G4), 


and  from  eqn  (G  28), 


= 


Cx  0  0  C: 

0  c-  C2  0 

o  c3  c4  o 

c3  0  0  Ci 


(39) 


*2  = 


/  0  0 

0  0  0 

GiC{lG:  -GtCi'CiG 3  /  0 

G;(C3-C4C;-'-Cx)G3  0  1 


(90) 


The  variables  used  in  eqns.  (39)  and  (40),  namely  C!C:iC3,C<Gt,C:G3  and  C4,  are 


defined  in  .\ppend;x  G. 
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Uses  of  the  Transfer  Matrices 

With  transfer  matrices  and  Xi  defined,  the  transfer  matrix  T  for  one  periodic  unit  cf 
the  planar  structure  can  be  obtained.  Figure  12  shews  the  planar  structure  which  has  been 
sectioned  into  four  basic  segments,  namely,  AlOOA'lAl ,  AKA'Kl'BlBLlAK,  BRB'„2'ClCL2BK  and 
CsCitlsCj,.  Ax A‘KYBLBL lAg  and  BkBI,2‘ClCi.2Bk  each  define  a  basic  periodic  unit  while 
AlOOA’lAl  defines  a  left  half-unit  and  C*C#33C*  defines  a  right  half-unit. 

Consider  first  section  AKA„  l'BLB'L\AK  in  Figure  12.  For  section  lt  —  l’L, 

[d, )  (d,  \ 


d,  Id, 

d"  ■  =Xl12  d"  ■  (41) 

Pn  Pn 

.P<  Jlt-lt  )ak-Ak 

where  X{n  is  the  transfer  matrix  which  includes  lengths  cf  i/2  in  both  main  members  I  and  II,  and 
where  \L  and  l‘L  are  points  just  to  the  left  of  points  1  and  1',  respectively.  Similarly,  for  section 
1*  ~1* .  where  1,  and  lj,  are  points  just  to  the  right  cf  points  1  and  1’,  respectively, 
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and  fa  section  3L  -  3'L , 
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.  du  ,  _  y 
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Substituting  eqns.  (41)  and  (42)  into  eqn.  (43)  gives 
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Since  section  AMA'K\BLBl\AK  is  representative  cf  a  periodic  unit,  the  transfer  matrix  T  for  one 
periodic  unit  is  obtained  from  eqn.  (44)  as 


(45) 


T  =Xln  XiXl11 

Similarly,  the  transfer  matrix  fcr  AL00'Ai_AL ,  denoted  by  T~u2  is  given  by 

T~ia  =  X}/2  X2  (46) 

and  the  transfer  matrix  for  C*C*3'3C* ,  denoted  by  Tv2,  is  given  by 

7 1/2  -X2X{/2  (47) 

Notice  that  the  transfer  matrices  7,  T~y2  and  7 1/2  are  all  transfer  matrices  which  relate  state 
vectors  from  left  to  right.  In  addition,  the  order  cf  the  transfer  matrices  on  the  right  hand  sides  cf 
eqns.  (45),  (46)  and  (47)  are  important. 

After  the  transfer  matrices  7,  7_1/2,  7 1/2  are  defined,  the  state  vectors  at  specific 
sections  of  the  planar  structure  can  be  obtained.  For  example  for  section  C  -C  , 
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As  another  example,  far  section  3-3', 
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The  frequency  response  functions  for  specific  sections  of  the  planar  structure  are  derived 
in  Appendix  H.  The  planar  structure  is  assumed  to  be  excited  by  sinusoidal  shear  forces  at  C  arxl 
C .  In  addition  to  deriving  the  frequency  response  functions  for  the  planar  structure,  Appendix  H 
also  gives  the  general  approach  whereby  the  problem  cf  forced  vibration  in  a  planar  structure  can  be 
treated.  From  eon.  (H9)  and  eqn.  (H10) , 
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where  the  superscripts  denote  the  excitation  location  and  the  subscripts  denote  the  response 
locations . 


From  eqn.  (H  11)  and  eqn.  (H  12), 

=  - 
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Fran  eqns.  (H  30),  (H  32)  and  (H  31), 
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\nriab!es  used  in  eqns  (50)  through  (56)  are  detuned  in  Appendix  H. 

The  impulse  response  functions  can  be  obtained  by  using  the  inverse  Fourier  Transform 

as  [10] 


W)  =3^  /  //;(«-)«'-</»  (57) 

where  the  subscript  q  denotes  the  response  location  a  interest  and  the  superscript  a  denotes  the 

excuiticn  location. 
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In  addition,  at  any  particular  frequency,  the  propagation  constants  can  be  found  by 
obtaining  the  negative  logarithms  cf  the  eigenvalues  cf  the  transfer  matrix  T  for  one  periodic  unit  of 
the  structure.  Notice  that  since  T  is  a  12x12  matrix,  there  are  six  pairs  cf  opposite  identical 
propagation  constants,  as  mentioned  earlier. 
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APPLICATION  OF  THE  TRANSFER  MATRIX  METHOD  TO  A  TETRAHEDRAL  TRUSS 

Figure  14  shows  a  tetrahedral  truss  which  is  used  in  this  section.  The  tetrahedral  truss 
is  assumed  to  consist  cf  identical  elastic  bars  with  distributed  mass  and  circular  cress  sections,  each 
of  length  (.  Each  connecting  bar  has  modulus  of  elasticity  E,  mass  density  p,  shear  modulus  G, 
cress -sectional  area  A ,  second  moment  cf  area  inertia  about  the  x  or  z  axis  J ,  second  moment  cf 
area  inertia  about  the  y  axis  JT  and  radius  cf  gyration  about  the  x  or  z  axis  i.  Furthermore,  the 
acrcss-sectional  dimensions  cf  each  bar  are  assumed  to  be  small  compared  with  its  length. 

State  Vectors 


For  each  connecting  bar  in  the  tetrahedral  truss  shewn  in  Figure  1,  the  state  vector  z 
consists  cf  six  force  components  and  six  displacement  components.  The  six  force  components  are 
M,M,T  y,y,  and  N ,  where  My  is  the  moment  about  the  y-axis,  M,  is  the  mem  ait  about  the  z- 
axis,  T  is  the  torque  about  the  x-axis,  V,  is  the  shear  force  along  the  z  direction,  V,  is  the  shear 
force  along  the  y  direction  and  N  is  the  axial  farce.  The  six  displacanent  components  include 
u  ,  v ,  w ,  dx,  »J»  and  0,  where  u  is  the  longitudinal  displacement,  v  is  the  lateral  displacanent  in  the  v 
direction,  w  is  the  lateral  displacement  in  the  z  direction  and  6,  v|i  and  8  are  rotations  cf  the  cress 
section  about  the  z,  y  and  x  axes,  respectively.  Thus,  for  each  connecting  bar, 


where  d  is  the  displacement  vector  such  that 


f:l 


* 


and  where  p  is  the  force  vector  such  that 


Figure  15  shows  the  sign  convention  for  the  force  and  displacement  variables. 

In  the  analysis  cf  the  tetrahedral  truss,  far  a  particular  station  (or  section),  the  state 
vectors  in  main  members  I  through  IV  are  needed.  The  reason  for  this  is  due  to  the  choice  of  the 
transfer  matrices  and  will  became  apparent  in  the  next  section.  Thus,  for  each  particular  station, 


where  the  subscripts  I,  II,  III  or  IV  denotes  the  member  number  for  which  d  or  p  is  defined  (refer 
to  Figure  14). 


Transfer  Matrices 


In  the  analysis  of  the  tetrahedral  truss,  three  transfer  matrices,  namely,  V,  V:  and  V3 
are  required.  Figure  16  shows  a  segment  of  the  tetrahedral  truss  and  Figure  17  shews  the  segment 
in  Figure  16  which  is  sectioned  along  planes  parallel  to  the  yz  plane  into  the  substructures 
responsible  for  transfer  matrices  V ,  V :  and  V3.  Referring  to  Figure  17,  the  sections  are  made  by 
curing  the  pericdic  unit  along  planes  GAkIDk,  ErK„FrLr  (cx  ElKlFlLl)  and  H3LJCL.  The 
subscripts  L  and  R  used  for  points  E,A,D,F,K,L,B  and  C  denote  points  just  left  and  right  of 
these  points,  respectively.  Thus,  the  first  transfer  matrix  V ;  represents  the  transfer  cf  state  vectors 
in  members  which  join  main  members  I  and  El.  For  example,  members  AD  and  BL  each 
constitutes  a  transfer  matrix  V';.  The  second  transfer  matrix  V’2  represents  the  transfer  cf  state 
vectors  in  four  bars,  each  of  length  t'2  in  main  members  I  through  IV  together  with  members 


which  connect  member  I  to  member  II,  member  HI  to  member  II,  member  I  to  member  IV, 
member  HI  to  member  IV.  Referring  to  Figure  17,  this  includes  members  ARKL ,  1FL  ,DrLl,  GEl  , 
AKfL,  DrFl  ,  A„El  and  DKEL .  The  third  transfer  matrix  V3  is  responsible  for  the  transfer  of  state 
vectors  in  four  bars,  each  of  length  (12  in  members  I  through  IV  together  with  member  which 
connect  member  II  to  member  I,  member  II  to  member  HI,  member  IV  to  member  I  and  Member 
IV  to  member  HI.  Referring  to  Figure  17,  this  includes  members  KkBl,  FkJ  ,  LkCl  ,  E„H  ,FMBL, 
E k  ,  E f f  B^  and  E ^ . 


The  transfer  matrices  VjVj  and  V3  are  derived  in  Appendix  I.  From  eqn.  (123),  the 
transfer  matrix  V  t  is  given  as 


V,  = 


0 


G1(Ci-CiC{lCl)G: 

0 

-GtC{lC  iGj 


0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

7  0  0  0  0 

u  u  0/000 

0  G£lC?Gl  0  0/00 

n  -  0  0  0  /0 

0  0  0  0  / 


0  0 
I  0 
0  / 
0  0 
0  0 


0  0 


0  C4G2-,G1 


Fran  eqn  (I  64),  the  transfer  matrix  V2  is  given  as 


(60) 


vs  = 


0 

c; 

0 

0 

0 

D, 


Z?3— C  3  D% 
D:+D,3^C}  d,  q 


0  0  0  0  0  Cj 

0  0  0  0  c2  o 

C3  0  0  c2  0  0 

0  c3  c2  0  0  0 

&:i  Ct  0  0  0 

£>6S0u-C3  D\2  0  C4  0  0 

O 

D: , 


0 

0 


0  c4  o 
0  0  ct 


161) 


Frcm  eqn.  (I  65),  the  transfer  matrix  V3  is  given  as 


0  0  0  0  c2 

o  0  0  Cj  0 

o  0  Cl  o  0 

C i  c2  0  0  0 

£i0+£I44C 3  c4  0  0  0 


£7+£i5-tC3  £13 


£2+£6+C3  £« 


|£j+£u+C3  £4 


0  C  't  0  0 

0  0  c4  0 
0  0  0  c; 


The  variable  used  in  the  transfer  matrices  in  eqns.  (60)  through  (62)  are  defined  in  Appendix  I. 


Uses  of  the  Transfer  Matrices 


With  transfer  matrices  V2iV 2  and  V3  defined,  the  transfer  matrix  T  for  one  periodic  unit 
of  the  tetrahedral  truss  can  be  obtained.  Figure  18  shows  the  tetrahedral  truss  sectioned  into  four 
basic  segments.  Referring  to  Figure  18,  segment  1  defines  a  left  half  periodic  unit ,  segment  2  and 
3  each  defines  a  periodic  unit  and  segment  4  defines  a  right  half  periodic  unit.  Following  the  same 
procedure  used  for  the  3-bay  planar  structure  in  the  last  chapter,  the  transfer  matrix  for  the  left 
half  periodic  unit ,  denoted  by  T  ~1/:,  is  given  as 

T<a  =V2V1  (63) 

Similarly,  the  transfer  matrix  for  a  periodic  unit,  denoted  by  T,  is  given  by 

T  =V2VtV3  (64) 

and  the  transfer  matrix  for  the  right  half  periodic  unit,  denoted  by  is  gisen  by 


=  V.V, 
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NUMER1CAL  EXAMPLE 

To  illustrate  the  application  of  the  transfer  matrix  method  in  the  analysis  cf  wave 
propagation  in  periodic  structures  the  case  cf  a  one-dimensional  elastic  rod  is  investigated.  The 
frequency  response  functions  for  specific  locations  for  the  longitudinal  vibration  in  an  elastic  rod  are 
obtained  Based  on  these  results,  the  impulse  response  function  for  a  location  in  the  rod  is 
generated,  whereby  the  wave  propagation  characteristics  far  both  square  pulse  excitation  and  a 
triangular  pulse  excitation  in  the  rod  are  studied. 

Figure  19  shows  the  rod  to  be  investigated  in  this  example.  The  rod  is  assumed  to 
consist  of  six  identical  rod  elements,  each  of  length  ■£,  with  no  impedance  mismatches  and  no 
material  damping  throughout.  The  material  and  geometric  properties  in  the  rod  are  given  as 


E  = 

7.46  X  10‘°Pa 

(lO&Xl&psi) 

A  = 

629  x  10-sm2 

(9.75x10  V 

P  = 

ZlClkg/m3 

(Oi  (bf/in3) 

(  = 

025  m 

(10  in) 

where  E  is  the  elastic  modulus,  A  is  the  cross -sectional  area  and  p  is  the  mass  density. 

The  frequency  response  functions  for  specific  locations  in  the  rod  due  to  a  sinusoidal 
axial  force  excitation  at  point  E  have  been  generated  in  eqns.  (231  through  (35).  Eased  on  these 
results,  a  basic  computer  program,  named  PROG1J3AS  (refer  to  Appendix  J),  is  written  to  obtain 
the  frequency  response  functions  numerically.  For  demons traticn  purposes,  a  plot  cf  the  frequency 
response  cf  the  force  at  B  due  to  excitation  at  E  versus  frequency  is  shown  in  Figure  20.  Jumps  in 
the  values  cf  trie  frequency  response  function  in  Figure  20  signify  resonance  conditions.  For 
comparison,  the  first  nine  natural  frequencies  for  the  rod  are  tabulated  in  Table  1  [14].  Modes  2.  5 
arc;  3  m  Table  1  are  net  shewn  m  Fig.  20  because  the  excitation  print  becomes  a  nodal  point  at 
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such  frequencies.  However,  since  PROG1JBAS  calculates  the  frequency  response  function  at  B  at 
discrete  frequency  intervals,  the  values  far  the  response  at  resonances  are  truncated  by  the  program 
and  do  not  reach  infinity.  The  impulse  force  response  function  at  B  due  to  a  unit  force  excitation  at 
E  is  generated  in  eqn.  (F44)  and  is  given  by 


MO  «  £  £  {-  *  -  (3+24  k)<y/jr]  -  -  (5+24  k)(\fjr) 

-ty  -  (7+24 k)€Vf]  -  <»  -  (9+24 

+  2(f  -  (15^-24*  +  <f  -  (17+24*  )*\/?] 

+  £(f  -  (19+24* +  S(r  -  (21+24*  )€*\/?]} 

E  E 


(66) 

Based  an  eqn.  (66),  a  plot  cf  the  impulse  farce  response  at  B  versus  elapsed  time 
intervals  is  shown  in  Fig.  21.  A  time  interval  of  which  is  numerically  equal  to  4.95  x  10'^ 

sec  is  used.  Note  that  t  Vf  is  the  time  required  for  the  impulse  excitation  to  travel  a  distance  i 

along  the  rod.  EXie  to  symmetry  conditions,  the  impulse  excitation  is  divided  into  equal  and  opposite 
going  waves,  each  having  magnitude  equal  to  half  a  Dirac  delta  function,  which  explains  the 
amplitude  of  the  response. 

Physically.  eqn.  (66)  states  that  if  a  unit  impulse  force  is  applied  at  E  at  time  zero,  a 
force  response  ct  amplitude  equal  to  half  a  Dirac  delta  function  will  be  observed  alter 

,  and  so  cn.  To  understand  this  further,  refer  to  F:g.  19.  At  time  zero,  a  unit  impulse 

force  cf  positive  magnitude  is  applied  at  E.  EXe  to  symmetry  conditions ,  this  force  is  divided  into 
two  equal  and  opposite  going  waves,  each  having  a  magnitude  equal  to  half  cf  the  applied 

amplitude.  After  three  time  steps  cf  the  left-going  wave  reaches  3.  Since  B  is  new  under 

compression,  it  experiences  a  force  cf  negative  magnitude.  This  explains  the  first  response  shown  in 


-3 


Fig.  21.  Now,  upon  reaching  3,  the  left  going-wave  trails  further  along  the  rod  until  it  reaches 


tte  boundarv  at  A  after  an  additional  The  wave  is  then  reflected  at  boundary  A. 

E 


However,  since  boundary  A  is  fixed,  there  is  no  sign  change  in  the  wave  due  to  reflection  and  after 
another  time  step  of  f\/^r,  it  reaches  B  again,  and  now  B  experiences  another  force  of  negative 

magnitude.  This  explains  the  response  at  in  Fig.  21.  Now  consider  the  right-going 


wave 


at  E.  After  two  time  steps  of  f\/-^  from  tirne  zero,  the  right-going  wave  reaches  boundary  G 

E  ^ 


and  is  reflected.  However,  since  boundary  G  is  a  free  boundary,  the  magnitude  of  the  wave  changes 
from  positive  to  negative.  After  another  five  time  steps  of  this  negative  wave  reaches  B . 


Now-  B  experiences  a  force  of  negative  magnitude.  This  explains  the  response  at  7f\/-£-  in  Fig. 

E 


21 .  The  negative  wave  then  travels  further  along  the  rod  and  after  it  reaches  boundary  A 


and  is  reflected.  Since  boundary  A  is  fixed,  there  is  no  change  in  the  sign  cl  the  wave  due  to 


reflection  and  after  another  <  V|  ,  it  reaches  B  and  B  again  experiences  a  negative  force.  This 


explains  the  response  at  in  Fig.  21.  Following  the  same  procedure,  each  individual 


response  at  3  as  given  by  eqn.  (66)  can  be  explained. 


Cnee  the  impulse  response  function  is  obtained  fer  a  particular  location,  the  response 
due  to  ether  forms  of  excitation  can  be  octained  using  the  relation  [10] 


.V'S(0  =  /  «AJ(f-T)  *  CO  d  (t)  (67) 

where  ,v*  (:;  is  tic  response  function  and  x(t)  is  the  excitation  function.  In  addition,  tor  the 
functions  (t)  and  ,h<  (r  — r),  the  subscript  a  denotes  the  type  of  response  (for  example,  axial 
force),  the  subscript  p  denotes  the  response  location  and  the  subscript  q  denotes  the  excitation 


■  ■  ' .  - .  .  -  . 
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location. 


For  the  present  case,  if  the  rod  is  excited  by  an  excitation  function  x(t),  the  response 
sy£  (t),  can  be  obtained  by  substituting  th ’  (r)  frcm  eqn.  (66)  into  eqn.  (67).  Because  (t)  in 
eqn.  (66)  consists  cf  eight  Dirac  delta  functions,  each  Dirac  delta  function  has  to  be  integrated 

separately.  Far  example,  consider  the  term  8 [t  —  (3+24*)-f\/-|-]  in  eqn.  (66).  The  response 

E 

v>’/  (0  due  to  this  function  is  given  by 


iy|(0  -  /  a[r-  (3+24*  )6  V?  “  t]  x  (r)dr 


(68) 


where  the  superscript  1  denotes  the  response  is  due  to  the  fist  Dirac  delta  function  in  eqn.  (66). 
Using  the  relation  [10] 

a 

/  8(r-r)/(T)rfT  »/(T)  (65) 

where  / (t)  is  any  function,  eqn.  (68)  becomes 

Ay/(0  **  [r-(3+24*)*VJ]  (70) 

The  same  procedure  can  be  applied  to  the  other  Dirac  delta  functions  in  eqn.  (66)  to  obtain  the 
total  response  sy§( t)  from  the  integral  in  eqn.  (67)  as 

.v>f(r)  =li{-*[:-3+24*)6V?l  -*[f'  ^24k)(^/jr] 

~  k  -o  ^  _ 

-x[r -<7+24*  )€%/?]  -x[t- (9*24Jfc)^V?l 
£  £. 

-r  [f -<15-24*  )l VJ}  -  x [r -(.17-24* )€ Vf  ] 

+r[r -(19+24*  )€\/£]  +  x[r-  (21+24*  )f V?]} 

(71) 

By  using  eqn.  (71),  the  response  at  E  can  be  obtained,  two  types  of  excitations  are 
considered  here  to  illustrate  the  approach. 

First  consider  the  rod  to  be  excited  at  £  by  a  square  raise  cf  unit  amplitude  and  cf 


>  l*-  .  ■  .  *  n  . 

-.■*  . .  - -- 


A 
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dura  cion  Thus,  the  excitation  Xi(r)  is  given  as 


*i(t)  = 


1  0  s  t  \/J 


ID 


otherwise 


(72) 


If  eqn.  (72)  is  substituted  into  eqn.  (71),  the  response  Ny§  (t)  is  obtained  as 


>yS( 0  -  ji {  (3+24*  )*Vf]  (5+24*)€\/f] 

(7+24*)€Vf]  -Xi[r-  (9+24*  )*V|-] 

+  X,[r-  (15+24* )*Vf]  +X,[r-  (17+24* )*VJ] 

+  X-.[r -(19.24*  )€V¥]  +X:[r-  (21+24*)*V^?]} 

£  E 


(73) 


Figure  22  shows  a  plot  cf  the  first  few  responses.  With  respect  to  time  intervals  of  i\/^r  each. 

E 

Similarly,  if  the  rod  is  now  excited  at  E  by  a  triangular  pulse  cf  duration  €\/-§-  and 

E 

peak  amplitude  cf  unity  such  that  X2(~)  is  given  as 


*:(t)  = 


T  V? 


•<fV? 


£  r  £ 
otherwise 


(74) 


4  -0 


If  sqn.  (74)  is  substituted  into  eqn.  (71)  the  response  v >1(0  is  given  as 

*>1(7)  =  (3-24*)fVf]  *Xj[/  -  (5-24*)tV|-] 

-Xj[r-(7+24*KVJ]  -X2[t-  (9+24*  )f\/J] 

+  X2[f -1(15+24*  K\^|]  +  Xj[f- (17+24*  K\/f] 
+  X:[f  -<19+24*  )<%/?]  +Xj[i-  (21J-24*)€\/?] 

E  E 


(75) 


n  :ure 


shows  a  plot  cf  the  first  few  responses  at  3,  with  respect  to  time  intervals  cf  i\'  ~ 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  transfer  matrix  method  is  observed  to  be  a  simple  and  straightforward  approach  in 
analyzing  wave  propagation  and  vibration  in  periodic  structures.  With  the  recent  advances  in  the 
field  of  transfer  matrices,  the  transfer  matrix  method  becomes  even  more  favorable.  The  transfer 
matrix  method  permits  a  simple  treatment  of  periodic  units  with  complicated  configurations. 
Moreover,  a  matrix  formation  is  most  suitable  for  periodic  structures  of  finite  total  length  since  the 
imposition  of  the  boundary  conditions  is  straightforward  as  is  observed  in  the  analyses  given  for  the 
one-dimensional  elastic  rod  and  the  3-bay  planar  lattice  structure.  Some  common  properties  cf 
transfer  matrices  are  outlined  in  Appendix  K,  which  may  serve  to  simply  calculations  in  applying 
the  transfer  matrix  method.  The  transfer  matrices  for  longitudinal  vibration  in  an  elastic  rod,  as 
given  in  eqns.  (B  11)  and  (B  16)  in  Appendix  B  and  the  transfer  matrix  for  flexural  vibration  in  a 
Timoshenko  Beam  as  given  in  eqn.  (D  16)  of  Appendix  D  are  used  in  Appendix  K  to  demonstrate 
the  properties  outlined. 

In  addition,  with  the  aid  cf  computers,  the  transfer  matrix  method  can  be  applied  with 
little  difficulty  once  the  transfer  matrix  for  a  periodic  unit  is  formulated,  as  is  demonstrated  in  the 
numerical  example. 

However,  there  is  one  setback  in  the  transfer  matrix  method.  If  a  periodic  structure  to 
be  analyzed  consists  of  a  large  number  cf  repeating  periodic  units,  the  numerical  application  of  the 
transfer  matrix  method  requires  multiplication  cf  a  large  number  cf  transfer  matrices  together,  the 
product  cf  which  may  become  enormously  large  and  difficult  to  handle,  even  with  the  aid  cf  a 
computer.  In  addition,  if  the  dements  In  a  transfer  matrix  are  frequency  dependent,  they  become 


larger  and  larger  with  increasing  frequencies.  The  accuracy  cf  the  transfer  matrix  method  may  be 
significantly  reduced  due  to  the  operations  cf  large  numbers.  This  is  demonstrated  by  the 
inaccuracies  in  obtaining  the  determinant  cf  the  transfer  matrix  for  flexural  vibration  in  a 
Timoshenko  Beam  at  high  frequencies  In  Appendix  X.  To  accommodate  for  this,  non- 


2- 


> 


C 


dime  nsi  ana  lized  forms  af  transfer  matrices  should  he  used.  For  illustration  purposes,  Appendix  L 
shews  the  derivation  cf  the  non-dimensionalized  farms  for  the  transfer  matrices  for  a  3-bay  planar 
lattice  structure.  Table  2  shows  the  calculated  values  cf  the  determinant  cf  the  transfer  matrix  cf  a 
Timoshenko  beam  (using  program  DET2.BAS  in  Appendix  K)  based  cm  various  values  cf  the  non- 
dimensionalized  parameter  flz  obtained  in  Appendix  L  using  both  single  precision  and  double 
precision  on  all  variables.  It  is  observed  that  the  value  cf  the  determinant  is  not  calculated  correctly 
to  be  unity  when  Cl2  reaches  a  value  cf  4  for  single  precision  calculations,  and  a  value  of  7 5  for 
double  precision  calculations.  Furthermore,  the  disagreement  in  the  calculated  results  far  single 
precisian  calculations  and  double  precision  calculations  again  justifies  the  inaccuracies  cf  the 
computed  results  due  to  operations  cf  large  numbers. 

Another  alternative  is  to  use  Cayley-Hamiltcn  theorem  [7]  when  multiplying  transfer 
matrices.  Cayley -Hamilton  theorem  expands  the  product  cf  a  number  cf  transfer  matrices  as  a 
linear  combination  cf  n  independent  analytical  functions  cf  the  transfer  matrix  T,  where  n  denotes 
the  dimension  of  T.  According  to  the  theorem, 


T‘ 


T‘~T~>  .  ,  T>  -T~> 

~Zal  — 5 -  +  bJ  - 5 - 

7-1  *  ^ 


(76) 


where  ajs  and  bjs  are  constants  and  K  is  the  number  of  times  transfer  matrix  T  is  multiplied  to 
itself.  The  constants  a, 's  and  bj 's  can  be  obtained  by  substituting  for  T  the  eigenvalues  cf  T  in  eqn. 
(72).  For  example,  if  transfer  matrix  T  is  a  4x4  matrix,  and  if  the  four  eigenvalues  cf  T  (refer  to 
Appendix  X)  are  substituted  successively  to  eqn.  (76),  four  independent  equations,  with  a/j  and 
bj's  as  the  only  unknowns,  are  obtained.  By  solving  these  four  simultaneous  equations,  the 
unknowns  ajs  and  bjs  can  be  determined.  Eqn.  (76)  can  even  be  applied  to  obtain  Tk .  The 
Cayley-Hamiltcn  theorem  becomes  favorable  when  the  number  of  multiplication  exceeds  the 
dimension  cf  the  transfer  matrix.  This  is  because,  according  to  eqn.  (76),  the  highest  power  in  T  is 
n,  2,  independent  on  the  value  a  k. 
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Table  2.  Computed  Values  cf  The  Determinant  cf  Transfer  matrix 
far  Timoshenko  Beam  far  Various  Values 
of  using  Double  and  Single  Precisions.* 


Determinant  of  Transfer  Matrix 
for  Timoshenko  Beam 


Elastic  modulus  is  7.46x  1010  Pt  (10.8X104  psi)  ,  shear  modulus  is  2.75x10'°  Pt  (4Dxl0*  psi), 
cross  sectional  area  is  6D4xl0_3m:  (9375 xl0_:  in),  second  moment  cf  area  inertia  is 
4 35x10"'°  in*  (1  £>9S  xlO'3  in*\  mass  density  is  7 2kg!m*(Q.l  Ibflin 3),  radius  cf  gyration  is 
2.7lx  IQ’3  m  Cl  J083SxlO~!  in),  arri  length  Ln  2.43x10-'  m  (9.75  in ) 


>  % 


Fig.  3  An  elastic  rod  loaded  by  a  sinusoidal  axial  -force  at 


TRANSFER  MATRICES 


MAIN  MEMBER 
H 


A  pianar  structure  sectioned  into  constituent  parts 
which  make  up  the  Transfer  matrices  Xj  and  X2. 


A  3-bay  planar  structure  sectioned  into  basic  periodic 
units  and  half  units. 


LENGTH  OF  EACH  BAR  *1 


MAIN  MEMBER  W 


Vyi  Myi  "  "  V 


-u,-*,T,N 


u,*,T,N 


Vy,  M  ,  f ,  v 


Vz »  0 , w 


Fig.  15  Sign  convention  for  the  forces  ond  displacements  in 
a  connecting  bar  of  a  tetrahedral  truss,  expressed 
in  local  coordinates. 
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LENGTH  OF  EACH  BAR  «l 


MAIN  MEMBER  ISO 

MAIN  MEMBER  HI 

MAIN  MEMBER  X 

MAIN  MEMBER  JL 


A  segment  of  the  tetrahedral  truss. 


transfer  transfer  transfer 

MATRIX  MATRIX  MATRIX 


TRANSFER 

MATRIX 


Fig.  17  Sectioning  of  the  segment  into  constituent  ports 
comprising  the  transfer  matrices. 
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Fig.  19  An  elastic  rod  clamped  at  the  left  end 
and  free  at  the  right  end,  loaded  by 
external  exitation  at  point  E. 
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APPENDIX  A 


PROPAGATION  CONSTANTS  AS  APPUED  TO 
THE  ANALYSIS  OF  WAVE  PROPAGATION  IN 
PERIODIC  STRUCTURES 


The  characteristics  of  wave  propagation  and  vibration  in  periodic  structures  are  best 
understood  in  terras  cf  propagating  and  non-prcpa gating  wave  motions.  In  genera],  waves  can  pro¬ 
pagate  in  some  frequency  bands  and  not  in  others  [2],  In  other  words,  periodic  units  in  periodic 
structures  behave  like  band-pass  filters,  responding  very  efficiently  in  certain  frequency  bands  only. 

Such  a  characteristic  is  commonly  described  by  a  propagation  constant  y  [12]  which  is  I 

described  by  the  nature  of  the  periodic  unit  and  the  corresponding  excitation  frequency.  The  har¬ 
monic  motion  at  one  end  cf  a  periodic  unit  is  equal  to  times  the  motion  at  the  other  end  from 
which  the  wave  is  travelling.  A  propagation  constant  y  can  be  real,  imaginary  cr  complex,  and  its 
value  always  occurs  in  positive  and  negative  pairs,  which  correspond  to  identical  but  opposite  going 
waves.  The  real  part  of  y  is  called  the  attenuation  at,  and  the  imaginary  part  is  called  the  phase 
constant  (or  wave  number)  k.  Purely  imaginary  propagation  constant  are  known  to  be  associated 
with  waves  which  propagate  energy,  whereas  purely  real  propagation  constants  belong  to  waves  of 
no  energy'  Cow  [1  ] .  The  frequency  bands  cf  the  real  part  cf  y  are  called  the  propagation  zones; 
aher  frequency  bands  axe  called  the  attenuation  zones.  The  number  of  possible  propagation  con-  | 

slants  (and  the  corresponding  waves),  for  a  periodic  unit  at  3  particular  frequency,  is  equal  to  twice 
the  number  of  state  vectors  (or  coupling  coordinates)  between  adjacent  periodic  units  [1],  For 
example,  for  the  flexural  vibrations  of  a  beam,  there  are  eight  propagation  constants  corresponding  | 

to  the  four  coupling  coordinates  which  are,  namely,  the  transverse  displacement,  the  angle  cf  rota¬ 
tion  of  the  cress  section,  the  shear  force  and  the  moment.  For  a  particular  value  cf  y  the  positive- 
going  waves,  each  as  a  function  cf  x  ,  and  of  the  fen  [2]  I 
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W,(x) 


2  a» 


(A 1) 


and  the  negative-going  waves,  each  as  a  function  of  x ,  are  cf  the  form  [2] 


W-(x)  =  2  Bne‘ 


(A2) 


Where  A,  and  B„  are  constants  and  IV  represents  the  wave  parameter  cf  interest  (for  example, 
stress  wave  in  a  periodic  structure). 


APPENDIX  B 


TRANSFER  MATRICES  FOR  LONGITUDINAL  VIBRATION 
IN  AN  ELASTIC  ROD  WITH  DISTRIBUTED  MASS 

In  this  appendix,  the  transfer  matrices  (from  left  to  right  and  from  right  to  left)  for 
longitudinal  vibration  in  an  elastic  rod  with  distributed  mass  are  derived  using  the  classical  wave 
equation  far  longitudinal  vibration  in  an  elastic  rod. 

Fig.  B1  shows  an  elastic  rod  with  distributed  mass  together  with  the  sign  convention 
adopted  far  the  farces  and  the  displacements.  The  rod  has  modulus  of  elasticity  E,  cross-sectional 
area  A ,  mass  density  p  and  length  €.  Furthermore,  N  (x  ,r)  represents  the  internal  axial  force  and 
u(x  ,t)  represents  the  longitudinal  displacement  in  the  rod. 

Consider  a  small  element  of  the  rod  as  shown  in  Fig.  B2.  Using  the  momentum  princi¬ 
pal, 


8.V 

dX 


-p 

dt1 


(Bl) 


By  the  definitions  of  stress  and  strain, 


N  =  uA  , 


3u  (x  ,r) 

t  =  ~i*  ’ 


and  c r  —  Et  , 


which  zive 


i£L  =AE^hilul 

dx  dx 2 


(B2) 


Equating  eqns.  (Bl)  and  (B2),  the  wave  equation  is  given  as 
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d2u(x,t)  E_  dlu  (x  ,t) 


dt 2 


dx1 


(B3) 


Now  let  u  ( x  ,t )  =  u  (x)sin(a*  +<{>)  where  w  is  the  circular  frequency  cf  vibration  and  cj> 
is  the  displacement.  Substituting  this  assumed  farm  into  eqn.  (B3)  gives 


-  a&inCa*  +$)  —  U  ^ - sin(ai  +4>)u  (x) 

dx*  E 


(B4) 


Cancelling  the  terms  — <j*sin(a*  -h*>)  from  eqn.  (B4)  gives 


iM£l+^e.u(j)=o 


CBS) 


Eqn.  (B5)  has  a  solution  of  the  form 

u{x)  - C isin(9y)  +C2  ccs  (0  j)  (B6) 

where  -'■V?  and  where  Ci  and  C2  are  constants. 

E 

transfer  matrix  from  left  to  right 

To  obtain  the  transfer  matrix  from  left  to  right  fa  a  rod,  adopt  the  left  end  cf  the  rod  in 
Fig.  B1  as  the  origin  for  the  x  axis.  Applying  the  boundary  conditions  to  eqn.  (B6), 

u  =  uL  ,ct  x  =  0 
u  =  ux  ,at  x  —  i 


where  uL  and  u,  are  known  quantities  and  solving  for  the  constants  in  eqn.  (B6) , 


u(x) 


ux  —jlcc$Q 
sin0 


sin(8j)  +  uiccs(0  j) 


(37) 


Using  eqn.  fB2)  cr  .V 


=  AE 


au 

dx 


and  calling  the  resulting  boundary  values  as, 


N  =Nl  ,  at  x  =  0 
and  N  =  Nk  ,atx=f, 


eqn.  (B7)  gives 


.  ,  i  sin9  ,, 

uK  =  cos  9  uL  +  —  -j-  Nl 


N/t  = - —  sin  9  Wi  +  cos  9  NL 


Writing  eqns.  (B8)  and  (B9)  in  matrix  form, 


tl- 


l  sin  9 
EA  9 


— jjbf  a?  008  ® 


where  ^  -  pA  and  the  subscripts  R  and  L  denote  the  right  and  left  ends  cf  the  bar,  respectively. 
Thus,  from  eqn.  (BIO) ,  the  transfer  matrix  T  from  left  to  right  is: 

f  i  sin  9 1 


cos  9 

EA 

e  2  sin  9 

9 

cos  9 

TRANSFER  MATRIX  FROM  RIGHT  TO  LEFT 

To  obtain  the  transfer  matrix  from  right  to  left  for  a  rod,  adept  the  right  end  cf  the  rod 
in  Fig.  B1  as  the  origin.  Then,  applying  the  boundary  conditions  to  eqn.  (B6), 

u  =  a*  ,  at  x  =0 
and  u  =  uL  ,  at  x  =■  —t  , 

where  u*  and  uL  are  known  quantities  and  solving  for  the  constants  in  eqn.  (36), 
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“00 


U*  COS  0— UL 

sin  6 


sin(8  j)  +  u,cos(8  y) 


Using  eqn.  (B2)  and  calling  the  boundary  values  as, 


eqn.  (B12)  gives 


N  =N’,  ,at  x  =  0 
and  N  =  NL  ,  at  x  =  —i  , 


uL  —  ccs  0  Uf  -  sin  8  SK 
FA  0 

Nl  =  sin  8  w,  +  cos  8  NK 


Writing  eqns.  (B13)  and  (B14)  in  matrix  farm, 


N 

v  J 


cos  8 


l  sin  8 
EA  8 


sin  8  cos  0 
0 


r  \ 
u 

N 


From  eqn.  (B15) ,  the  transfer  matrix  T~l  from  right  to  left  is 


_  J_ 

cos  0  EA 

s*n  9  ros  9 


sin  8 
8 
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APPENDIX  c 

TRANSFER  MATRIX  FOR  LONGITUDINAL  VIBRATION  IN  AN  ELASTIC  ROD 
WITH  DISTRIBUTED  MASS  AND  DAMPING 


In  this  appendix,  the  transfer  matrix  from  left  to  right  for  longitudinal  vibration  in  an 
elastic  rod  with  distributed  mass  and  damping  is  derived  using  the  wave  equation.  The  derivation  d 
the  transfer  matrix  from  right  to  left  in  an  elastic  rod  with  distributed  mass  and  damping  follows  a 
similar  approach  and  therefore  is  not  discussed. 

Fig.  Cl  shows  an  elastic  rod  with  distributed  mass  and  damping  together  with  the  sign 
convention  adopted  far  th:  forces  and  the  displacements.  The  rod  has  modulus  ci  elasticity  £, 
material  damping  c,  mass  density  p,  a  css -sectional  area  A  and  length  €.  Furthermore,  N(x  ,t) 
represents  tte  internal  axial  force  in  the  rod  and  u  (*  ,r )  represents  the  longitudinal  displacement  in 
the  rod. 

Consider  a  small  element  of  the  rod  as  shown  in  Fig.  C2.  Using  the  momentum  prinri- 
ple, 


New  assume  that  tie  rod  can  be  modeled  as  a  simple  Vcigt  material  [11]  with  elastic  modulus  E 
and  damping  constant  c  such  that  the  stress-strain  relation  is  given  by 


(C2) 


3v  the  definitions  cf  stress  and  strain. 


since  N  =  a A 


and  e 


3u  (x  ,t) 
dx 


dN 

dx 


=  A 


r.  d2u  (x  . t  ) 
E~ ~ 


dl 

.  t*1  . 


(C3; 


Equating  eqns.  (Cl)  and  (C3)  and  rearranging, 


dzu  (x  ,f  ) 

.  =1 

r.d2u(x,t)  |_a  [a«J(x,r) 

\  ' 

dt2 

P 

dx2  dt  dx2 

(C4; 


Now  assume  that 


u(x  ,t )  =u(x)e/t 


(cs; 


Alter  substituting  eqn.  (C5)  into  eqn.  (C4)  and  cancelling  the  e‘m  terms, 


M(x)  =0 

ax2  (£+ic)  v  ; 


(cs; 


Eqn.  (C6)  is  the  wave  equation  for  longitudinal  vibration  in  an  elastic  rod  with  distri¬ 


buted  mass  and  damping.  It  has  a  solution  cf  the  farm. 
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«0)  -Ci  sin  (6j)  +  C2  cos  (0y) 


(C7) 


where  0  =  •fai'v  — 
v  E  +ic 


and  there  C t  and  Cj  are  constants. 


Applying  the  boundary  conditions, 


u  =  uL  ,  at  x  =0 
and  u  =uK  ,  at  x  =  ( 

where  uL  and  un  are  known  quantities  and  calling  the  axial  forces  at  the  left  end  and  the  right  end 
of  the  rod  SL  and  .V, ,  respectively,  from  eqn.  (C7), 

U„  =  cos  fc,  +  (CS) 

FA  ft 

and  N„  = - sin0ut  +  cos  WL  (C9) 

Writing  eqns.  (C8)  and  (C9)  in  matrix  form, 


.  (  sin0 

u 

0X9  EA  8 

u 

V 

X 

COS  0 

N 

K  > 

(CIO) 


where  u  =  pA . 
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APPENDIX  D 

TRANSFER  MATRIX  FOR  FLEXURAL  VIBRATION  IN  AN 
ELASTIC  BAR  INCLUDING  THE  EFFECT  OF  SHEAR  DEFLECTION 
AND  ROTARY  INERTIA 


The  transfer  matrix  for  flexural  vibration  in  an  elastic  bar  including  the  effect  of  shear 
defection  and  rotary  inertia  is  derived  in  this  appendix.  Fig.  D1  shows  an  elastic  bar  with 
distributed  mass  and  the  sign  convention  chcsen.  The  bar  has  modulus  of  elasticity  E,  shear 
modulus  G ,  second  moment  cf  area  inertia  about  the  y  axis  J ,  mass  per  unit  length  u,  and  radius  of 
equation  i.  In  addition,  w  denotes  the  lateral  displacement,  <j>  denotes  the  rotation  cf  the  crc&s- 
sectional  area  about  the  y  axis  and  M  and  V  represent  the  moment  and  the  shear  force, 
respectively. 

Consider  first  an  element  of  the  bar  as  shown  in  Fig.  D2,  which  gives 


V 


(Dl) 


where  GAs  =  G AIK,  is  the  shear  stiffness  and  K,  is  the  form  factor  which  depends  on  the  shape  of 
the  crcss-sectional  area. 


The  constitutive  bending  relation  for  a  bar  is 


M  =EJ  4^ 
ax 


(D2) 


Now  consider  Fig.  D3,  which  shows  a  free-body  diagram  of  an  element  cf  the  bar. 
Equilibrium  consideration  give  the  following  equations: 


V  -  (D3) 

-  (D4) 

Differentiating  eqn.  (Dl)  with  respect  to  x,  using  the  relationship  in  eqn.  (D2)  and  substituting  into 
eqn.  (D4)  gi\es 


dx 

dV_ 

dx 


dhv 

dx2 


+ 


u tu 

GA, 


w 


M 

EJ 


0 


(D5) 


Now  differentiating  eqn.  (D3)  with  respect  to  x ,  using  the  relationship  from  eqn.  (D2)  and  then 
substituting  into  eqn.  (D4) , 


dx2  EJ  ^ 


Eliminating  M  in  eqns.  (D5)  and  (D6) 


d*w  crHr  d2w  g*-QT 

dx2  e  dx2  ~  v 


,  PIlJ2^2 

where  a  =  — 

GA, 


(D6) 


(D7) 


EJ 


and  3* 


jjfaf  t'4 

EJ 


Now  assume  a  displacement  w  such  that 


w  =CeXjlt 


where  C  is  a  constant.  Substituting  the  assumed  form  far  displacement  into  eqn.  (E7),  cancelling 
the  term  -CeKllt  from  the  equation,  the  characteristic  equation  is  obtained  as 


X4  +  (cH-t)X2  —  ({3* -err)  =  0 


The  roots  cf  eqn.  (D8)  are  ±  ^  and  ±  j'Xj  where 


V?  +  J  (°^0  ?  J  (w-r) 


with  X|  —  Xi  =  <r  +  t 


X?X2  =  0*  -  err 


Therefore,  the  solution  is 


w  +  C2e-M"‘  +  C3eJ*I,'t  +  C4e 


j.  r  J**'* 


where  Ci,  C2,  C3  and  C4  are  constants. Eqn.  (D9)  can  be  written  in  the  form 


w  =CtCCsh(Xi  j)  +C2sinh(^  j) 
+  C3cos(X2  j)  +  C4sin(X2  j) 


(DIO) 


where  Ci  =  C i  +  C2 


C2  —  Ci  —  C2 


C3  -  C3  +  C4 


C4  =j(C3~Ct) 


Examination  cf  eqns.  (Dl)  and  (DIO)  shows  that  V  and  w  are  of  the  same  form. 


Therefore,  let 


V  —  /41cosh(X1  y)  +-42Sinh(Xi  y)  +  M:jCCs(X2  y)  +  ^in^ 


Using  eqns.  (Dl),  (DIO)  and  (Dll), 


* 


e 

&EJ 


|(<r+\i)[A1cosh(X1  y)  +  A2sinh(X1 


+  (cr-X|)  [A 3ccs(X2  y)  +  A4sin(X2  y)] 


Using  eqns.  (D2),  (DIO)  and  (Dll), 


(cr-*-X?)  -y  [^sinh^  y)  +  /42ccsh(^ 


-  (<r-X|)  y-  [Ajsmfo  ^) 


' 


Writing  eqns.  (DIO)  to  (D13)  in  matrix  form, 


/ 


€J(cr+X?)  x 

— — - ccsh(X1  -r 

$EJ  t 

i\{c H-Xp)  _._L/X  X, 
- ^ - sinh(Xl  y) 


ccshfXi  ”) 


*X, 


i,  (cH-X?) 


(3* 

sinh(X! 
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~j£E7  f >  |eT  “<>*  f ! 

cosfk  i)  ilfctll  si„(>s  i) 
p^y  v  2  t  *  $ej  K  i  * 

€2\i((r-\i)  x  s  ^^2(cP“X|)  x 

- pT - softy)  - ^ - cds(X21 

®s(X2  j)  sin(Xi  j) 


(D 14) 


Imposing  the  boundary  values, 


w  =  wL  ;  «Ji  =  4»t  ,  at  x  =  0 


w  =wM  ;  ij/  =  4/j,  ,atx~i  , 


wtere  wL ,  w„ ,  ,  ifo ,  are  known  quantities.  Calling  the  shear  forces  and  moments  at  the  right  end 

and  the  left  end  cf  the  beam  as  V,  and  MK,VL  and  ML ,  respectively,  eliminating  the  constants 
A2,  A3,  and  At,  the  transfer  matrix  is  obtained  which  relates  the  state  vectors  at  x  -  0  and  x  =  l 


Cq—VC  2 

rCi 

&EJ 


<[c  1  -  (o—t)c3] 

Co  -  TCj 

-y-  [-TC!  +  (3‘-'-T2)c3] 
:j)  cl  Cl 


m 


s 


[“*1  +  (P*+oz)c3] 


?c2  e 


EJ  3 *EJ 

t(Cl~TCy)  V- 

EJ  EJ  Cl 
c0  -  rC2  {[ ci  ~  (cHr)e3] 


c0  -  cr2 


with  A  =  — - - 

\2 


c0  —  i\(X|ccBhXl  +  X^cqbXj) 


Cx  =  GXhh  +  T~  sin As) 

Ax  Aj 


c2  =  A(cosh\x  -  cosXj) 


c3  »  A( 


sinAXx  sinXj 


) 


From  Equation  (D15),  the  transfer  matrix  T  required  is  thus 


Co  -  CC; 

3* 

TC3 

&EJ 

Jhrci 


t[c I  -  (o^)Cj] 
Co  -  TCj 
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£c_2 

EJ 

t(cx 


rca) 


EJ 

Co  -  tc2 


€ 


Cj 


P’EJ' 

EJ  Cl 

€[ct  -  (<M-r)cj] 
Co  ~  «2 


(D16) 


z 


Fig.  D1  An  elastic  bar  with  distributed  mass  (Timoshenko 
beam  model). 
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APPENDIX  E 


TRANSFER  MATRIX  FOR  LONGITUDINAL  AND  FLEXURAL  VIBRATION  IN 
AN  ELASTIC  BAR  WITH  DISTRIBUTED  MASS  AND  ROTARY  INERTIA 


The  transfer  matrix  for  longitudinal  and  flexural  vibration  in  an  elastic  bar  with  distri¬ 
buted  mass  and  rotary  inertia  is  derived  in  this  appendix.  Fig.  El  shews  an  elastic  bar  with  the  sign 
convention  adopted  for  the  forces  and  displacements  in  the  bar.  The  bar  has  modulus  of  elasticity 
E ,  shear  modulus  G ,  second  moment  of  area  inertia  about  the  y  axis  J,  radius  cf  gyration  about  the 
y  axis  i,  crcss-sectianal  area  A  and  length  t.  Furthermore,  the  bar  has  displacement  components 
which  consist  of  transverse  displacement  w,  longitudinal  displacement  u,  rotation  cf  the  bar’s  cross- 
section  \J»  and  farce  components  which  include  the  shear  force  V ,  axial  force  N  and  moment  M . 

If  the  transverse  deflections  are  assumed  small  (relative  to  the  bar  cross-section),  the 
longitudinal  and  flexural  vibrations  (or  waves)  are  not  coupled.  The  transfer  matrix  can  then  be 
obtained  by  directly  assembling  the  transfer  matrix  fer  longitudinal  vibration  with  the  transfer 
matrix  for  flexural  vibration. 

Since  for  longitudinal  vibration  in  an  elastic  bar  with  distributed  mass, 
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and  for  flexural  vibration  cf  an  elastic  bar  with  distributed  mass  and  rotary  inertia, 
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where 
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Thus,  the  transfer  matrix  T  specified  is 
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APPENDIX  F 


FREQUENCY  RESPONSE  AND  IMPULSE  RESPONSE  FUNCTIONS  FOR 
LONGITUDINAL  VIBRATION  IN  AN  ELASTIC  ROD 

Figure  FI  show*  the  elastic  rod  to  be  used  in  this  Appendix.  The  rod  has  modulus  cf 
elasticity  £ ,  mass  density  p  and  crcss-sectianal  area  A .  N  denotes  the  axial  force  and  u  denotes  the 
longitudinal  displacement  in  the  rod.  For  demonstration  purposes,  the  rod  is  damped  at  one  end 
and  the  other  end  is  left  free.  The  rod  is  loaded  at  Section  E  with  a  sinusoidal  axial  farce  cf 
magnitude  N0 .  Furthermore,  the  rod  is  assumed  to  consist  cf  six  identical  rod  segments ,  each  of 
length  (.  In  this  Appendix,  the  frequency  response  functions  at  sections  A,  B  ,C ,  D ,  E  and  F  will 
be  obtained  using  transfer  matrices.  After  the  frequency  response  functions  are  generated,  the 
impulse  response  functions  can  be  obtained  using  inverse  Fourier  transforms.  In  particular,  the 
impulse  response  function  is  obtained  for  section  B. 

FREQUENCY  RESPONSE  FUNCTIONS 

For  longitudinal  vibration  in  an  elastic  rod,  the  state  vector  r  consists  cf  the  longitudinal 
displacement  u  and  the  internal  force  N .  Thus, 


At  boundary  A, since  the  displacement  is  defined  (that  is,  u  =  0), 
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Ai  boundary  G ,  since  the  force  is  defined  (that  is,  S  =  0), 


Tie  excitation  is  given  by  Le‘~  where  I  is  an  input  vector  whose  dements  correspond  to 
those  in  the  state  vector  being  considered.  Thus,  in  general. 


L  = 


(F4) 


where  D  is  the  displacement  excitation  and  P  is  the  farce  excitation.  In  the  present  case,  since  the 
axial  force  at  £  is  the  excitation, 
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For  longitudinal  vibration  in  an  das  tic  rod,  the  transfer  matrix  relates  the  state  vectors 
at  the  two  ends  of  a  rod  segment  as,  from  left  to  right, 
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and  from  right  to  left, 
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Let  the  transfer  matrices  giren  in  eqn.  (F6)  and  (F7)  be  T(()  and  T'1^),  respectively, 
wtere  the  t  in  parentheses  signifies  that  the  transfer  matrices  are  for  a  rod  cf  finite  length  (. 
Notice  that  since  there  is  no  impedance  mismatch  between  the  individual  sections  and  the  rod  can 
be  considered  as  one  continuous  member,  multiplying  T  (£)  rt  times  is  equivalent  to  replacing  €  in  T 
by  rti,  or 


T-(<)  =T(n() 


For  the  state  vector  at  G , 


{o}G  =r(f)r(<)r(f)r(f)T(f)r(f)  +  r(«)r(€)  . 


or,  using  eqn.  (FS),  eqn.  (F9)  can  be  written  as 
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Now  let  A  =  T (6f)  where  A  is  a  2 >2  matrix  with  elements  an,  al2,a2i  and  au  and  let 
B  =T (2f)  where  B  is  a  2X2  matrix  with  elements  bn,  bn,  b2 \  and  bn-  Eqn.  (F10)  can  be  written 
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Considering  the  force  vectors  only,  from  eqn.  (Fll), 
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from  which 
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examination  ctf  eqn.  (F6)  gives 


an  =  cos69 


and  bn  —  cos29 


Substituting  eqns.  (F13)  and  (F14)  Into  eqn.  (F12), 
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Sinoe  by  definition  H  (w) 
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where  the  superscript  E  denotes  the  exitation  location,  the  subscript  A  denotes  the  response  location 
and  the  subscript  N  denotes  a  force  response .  .At  g 


(FI  7) 


Let  c  .i  ,  d2,  c:i .  c::  be  elernents  cf  the  2  >2  matrix  7  (£)  such  that  eqn.  (F17)  becomes 
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Using  eqns.  (F15)  and  (FIS), 
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Examination  rf  eqn.  (F6)  gi\es 
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from  which, 
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Notice  that  when  there  is  no  subscript  to  the  left  cf  the  frequency  response  function  is  a 

vector  whose  top  element  represents  the  displacement  response  and  whose  bottom  element  represents 
the  force  response. 

Similarly,  for  station  C  , 
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for  station  D, 
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and  for  station  EL  where  EL  is  the  point  just  to  the  left  cf  E , 
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Now  the  frequency  response  functions  at  G,  F  and  Ex  are  obtained  where  £>  is  the 
point  just  to  the  right  of  E.  To  avoid  crossing  the  discontinuity  due  to  the  excitation  at  E ,  transfer 
matrices  from  right  to  left  are  used. 

Using  Equation  (F7) ,  for  the  state  vector  at  A , 
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The  minus  sign  is  due  to  the  sign  convention  chosen.  This  is  because  transfer  matrices  from  right  to 
left  is  used.  Now  let  D  =  T~l(6()  with  elements  dn,  d d2i,  and  dn  and  E  =7*-1(4f)  with 
dements  en,  e;:,  e21  and  e-^.  Equation  (F26)  can  be  rewritten  as 
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Considering  only  the  displacement  vector  at  A , 
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From  eqn.  (F7) 
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Thus, 
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where  the  subscript  u  denotes  displacement  response.  At  station  F , 
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Using  eqns.  (F28),  (F29),  (F30)  and  (F33), 
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Since  eqn.  (F7)  gi\es 
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eqn.  (F34)  can  be  written  as 
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Similarly,  for  station  £*  where  £,  is  just  to  the  right  cf  point  £, 
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IMPULSE  RESPONSE  FUNCTIONS 

In  the  present  case,  since  the  frequency  response  functions  are  quite  simple  (from  a 
mathematical  point  cf  view),  the  impulse  response  functions  can  be  obtained  in  dosed  form  by 
simply  taxing  the  inverse  Fourier  transform  cf  the  frequency  response  functions. 

First  consider  the  impulse  force  response  functions  at  point  3.  Frcm  Equation  (F22), 
for  the  force  response, 
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where  the  subscript  ,V  denotes  the  force  response.  Using  the  relationship 


Equation  (F37)  can  be  rewritten  as 
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Rearranging  eqn.  (F38), 
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Multiplying  eqn.  (F39)  by  - —  and  arranging  terms, 
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Factonng  the  term  from  eqn.  (F40) 


,-::a 


co)  =  ~-{e‘a' 


i 


Since  - - - - ^  e‘:4*  [13],  eqn.  (F4 1 )  becomes 
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Rearranging  terms  in  eqn.  (F42) 
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Since  MO  =  ^  /  //(u»)«'-<icu  [10],  and  —  /  e^e-du 
the  Dime  delta  function,  also  keeping  in  mind  that 
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6(r  — r)  where  £(r  — r)  is 


integration  cf  eqn.  (F43)  term  by  term  gives 
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APPENDIX  G 

TRANSFER  MATRICES  FOR  WAVE  PROPAGATION  IN  A  3- BAY 
PLANAR  LATTICE  STRUCTURE 

Fig.  G1  shew  the  3-bay  planar  lattice  structure  to  be  used  in  this  appendix.  The  planar 
structure  is  assumed  to  consist  of  identical  elastic  bars  with  distritxited  mass,  each  of  length  l. 
Each  bar  has  modulus  of  elasticity  £ ,  mass  density  p,  shear  modulus  G ;  cross-sectional  area  A , 
second  moment  cf  area  inertia  about  the  y  axis  J  and  radius  of  gyration  about  the  y  axis  i.  The 
cress -sectional  dimension  cf  each  bar  is  assumed  to  be  small  compared  to  the  length,  and  the 
Timoshenko  beam  model  is  adopted.  Each  bar  has  transverse  displacement  w,  longitudinal  displace¬ 
ment  u ,  angle  of  rotation  of  the  cross-section  <Jr,  axial  force  N ,  shear  force  V  and  taanent  M .  The 
sign  convention  for  the  displacements  and  forces  in  a  connecting  bar  as  shown  in  Fig.  G2. 

In  analyzing  the  wave  propagation  and  vibration  cf  the  planar  structure  shown  in  Fig. 
Gl,  two  transfer  matrices  are  involved.  The  first  transfer  matrix  X2  involves  the  transfer  cf  state 
vectors  in  two  bars,  each  cf  length  t  in  main  members  I  and  n.  For  example,  members  12  and  12’ 
together  are  represented  by  such  a  matrix.  The  second  transfer  matrix  X2  involves  the  transfer  cf 
state  vectors  across  the  junctions.  In  a  periodic  unit,  the  members  that  join  main  members  I  and  II 
constitute  such  a  transfer  matrix.  Far  example,  member  IT  results  in  a  transfer  matrix  X2.  Fig. 
G3  shews  a  3-bay  planar  structure  which  has  been  sectioned  into  its  constitutive  parts  responsible 
for  transfer  matrices  Xt  and  X2.  The  subscripts  R  and  L  are  used  to  denote  points  which  are  just 
to  the  right  and  points  which  are  just  to  the  left,  respectively,  cf  junctions  which  connect  main 


members  I  and  U. 
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transfer  MATRIX  Xi 

Transfer  matrix  can  be  obtained  by  simply  assembling  the  transfer  matrices  for 
flexural  and  longitudinal  vibrations  in  straight  ban  such  that  both  the  state  vectors  in  main 
members  I  and  II  are  considered.  The  transfer  matrix  for  longitudinal  and  flexural  vibrations  in 
straight  ban  is  given  as 
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Partitioning  the  transfer  matrix  given  into  four  ^x3  submatrices  CUC 


such  that 
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2)C3,  and  C, 


The  following  relations  can  be  written 


where  the  subscripts  I  and  II  stand  for  member  numbers. 


Eqns.  (Gl)  and  (G2)  can  be  combined  to  form 
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Notice  that  eqn.  (G3)  relates  the  state  vectors  at  the  right  end  to  that  cf  the  left  end  cf 
a  section  of  a  periodic  unit  which  is  represented  by  transfer  matrix  Xx. 

From  eqn.  (G3), 


TRANSFER  MATRIX  X2 
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From  the  previous  section, 
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Consider  an  I-juncticm  as  shown  in  Fig.  G4.  The  forces  and  displacements  cf  member  11'  (Fig.  G4) 
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are  shown  in  Fig.  G5.  For  member  11', 
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where  d  and  p  are  displacement  and  force  vectors  in  the  local  11'  coordinate. 

In  crcx sing  over  from  left  to  right  of  junctions  1  and  1',  the  displacements  are  unchanged 
but  the  bar  11'  introduces  discontinuities  in  the  forces  and  moments.  Assuming  the  displacements  at 
the  junctions  are  known,  these  forces  may  be  computed  from  the  elastic,  geometric  and  mass  proper¬ 
ties  of  bar  11  ’. 

Fig.  G6  and  G7  show  the  directional  relationships  between  local  and  global  state  vectors 
in  junctions  1  and  1',  respectively.  By  inspection  cf  Fig.  G6, 
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Similarly,  inspection  cf  Fig.  G7  gives 
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From  eqn.  (G5), 
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Multiplying  eqn.  (G14)  by  C{1  and  rearranging, 
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Transforming  the  local  coordinates  to  the  global  coordinates  in  the  main  members,  using  eqns 
(G10)  to  (G13), 
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Writing  eqns.  (G17)  and  (G18)  in  matrix  form, 


fc)- 


G2C4Cj-1G1 

GtC{lGi 


G2(C3  -c4c2-1c1)g, 

-G4C2_1C1G3 


(G19) 


New  examine  Eg.  G8,  which  shows  the  forces  at  junction  1.  From  equilibrium  con¬ 


siderations, 
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where  the  subscripts  L  and  R  stand  for  the  left  and  right  end  points  of  1.  Similarly,  for  junction  1', 
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From  eqns.  (G20)  and  (G21), 
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where  /  is  the  identity  matrix,  and  from  eqn.  (G19), 
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As  mentioned  previously ,  displacements  are  continuous  in  crossing  junctions  1  and  1',  thus, 
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'  \ 

/  0  0  0 

d\ 

0/00 

dv 

U-j 

l = 

0  0  0  0 

Pv 

0  0  0  0 

P\ 

\ 

i  ,  4 

(G23) 


Eqns.  (G23)  through  (G25)  can  be  combined  to  fcnn 


Pv  GjCiGi  -G.Cj-'CiCj  1  0  Pv 

P\\K  G2C4C2~lG1  G2(C2l-C4Cf1G1)G3  0  I  (pi 


Keeping  in  mind  that  1  ccrres ponds  to  a  point  in  member  I  and  1  ’  corresponds  to  a  point 
in  member  n,  eqn.  (G26)  can  be  rewritten  as 


/  0  0  oM' 

0  I  0  0  d„ 

G4C{xGx  -C*CilCiC3  I  0  Pn 

G£4CilG2  G2{C2l-C4CixGx)G3  0  1  {P> 


(G27) 


Notice  that  eqn.  (G27)  relates  the  state  vectors  (on  main  members  I  and  II)  an  the  right 
end  to  the  state  vectors  an  the  left  end  of  a  section  af  a  periodic  unit  which  is  represented  by 
transfer  matrix  X2. 

From  eqn.  (G27), 


g4C{xg2  -gac;'c2c  3  /  o 

CrC4Cj-lG:  G2(C2.-C 4C{lG:)G3  0  I 


(G2S) 
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Fig.  G3  A  planar  structure  sectioned  into  constitutive 
parts  which  make  up  the  transfer  matrices 
X1  and  Xg. 


Directional  relations  between  global  and 
local  state  vectors  at  location  1. 
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Directional  relations  between  global  and 
local  state  vectors  at  location  V. 
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Fig.  G3  Forces  at  location  1 
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APPENDIX  H 

FREQUENCY  RESPONSE  FUNCTIONS  FOR 
A  3-BAY  PLANAR  LATTICE  STRUCTURE 

The  frequency  response  functions  for  specific  locations  in  a  planar  structure  can  be 
obtained  through  the  use  of  transfer  matrices. 

Consider  the  3-bay  planar  structure  shown  in  Fig.  HI.  For  demonstration  purposes,  the 
structure  is  loaded  at  the  mid  points  between  members  23  and  2 ’3’  with  sinusoidal  shear  forces.  It 
will  be  shown  that  the  frequency  response  functions  for  locations  A ,  A',  B,  B',  4,  5  and  6  can  be 
ob faired.  For  simplicity,  assume  that  the  crcss-secticnal  dimensions  in  the  connecting  bars  are 
small  compared  to  tte  lengths  and  that  the  structure  consists  cf  identical  rod  elements  throughout. 

The  excitation  is  given  by  Le1*  where  L  is  an  input  vector  whose  dements  correspond  to 
those  in  the  state  vector  being  considered.  In  ether  words, 


where  D  is  the  displacement  excitation  vector  which  includes  the  longitudinal  displacement  u,  lateral 
displacement  w  and  rotation  cf  the  bar  crcss-secticnal  area  i>,  and  ?  is  the  force  excitation  vector 
which  indudes  the  moment  M ,  shear  force  V  and  amal  feres  N.  The  subscripts  I  and  II  used  in 
eqn.  (HI)  stand  for  the  entire  upper  and  lower  member  numbers.  In  the  present  case,  since  the 
shear  forces  are  the  only  excitations, 


and  at  3-3 


Now  let  be  the  transfer  matrix  corresponding  to  two  bar  elements  cf  length  (  in 
members  I  and  II,  Xm  be  the  transfer  matrix  corresponding  to  two  bar  elements  of  length  €/2  in 
members  I  arri  n,  arei  X2  be  the  transfer  matrix  fcr  the  vertical  members  which  connect  main 
members  I  and  D.  For  the  state  vector  at  section  3-3’,  going  from  left  to  right  along  the  structure, 


0  '  =  x2x1x2xlx2xlx2 0  +x2x}/2-L2 

0  0  i, 


Let  K  =  X:X,X:X  X:X  X:  be  a  12x12  matrix  with  elements  cf  3x3  submatrices  where  each 
submatrix  is  denoted  by  k,t  (i  =  12.3,4  and  j  —  1234)  and  let  B  —  X2X\!2  be  a  12x12  matrix 
with  3x3  submatrices  b,,.  Then  considering  only  the  lead  vectors  in  3-3’,  from  eqn.  (H4), 


At  station  A -A', 
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*31*22 

*41*42 


*33*14 

*43*44 


Li)c  -C' 


(H 5) 


> 

d, 

d, 

dji 

N 

* 

II 

djj 

P, 

0 

Pi 

0 

A  -A' 

4  J 

Let  B'  =  XlnX i  WiJi  elements  cf  3x3  submatrices  b,),  then 


'  N 

d, 

*11*12 

d„ 

*2i*n 

Pn 

b  31*32 

Pi 

J 

A  -A ' 

*41*42 

d‘\ 

&U  Jo-O' 


Substituting  eqn.  (H5)  into  eqn.  (H6), 


(H6) 


4 

bn 

*u 

bn 

*22 

P// 

b  3i 

bn 

/>/ 

/ 

A  -A' 

*« 

b’,2 

■ 

-1 

■ 

*31*32 

*33*34 

^2 

*42*42 

*43*44 

Li 

• 

-  / 

Since  H  (m)  is  bv  definition  -resP?*se~-  from  eqn.  (K7), 

excitation  ^  v 


(H7) 


*u 

*22  1 

.1 

■ 

*21 

*n 

*31*33 

*33*34 

*31 

*n 

*41*42 

*43*44 

*i. 

*i2 

• 

where  the  subscript  in  H  (<a)  stands  for  the  response  location  and  the  superscript  stands  for  the 
excitation  location.  Sectioning  eqn.  (HS). 
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*«1  *42  J  1*41  *42  I  1*43  *44 


bn  b  i 


=  ~ 


*21  bzi  *31  *32  *33  *34 

*31  *32  *  41  *42  *  43  *  44 


(H10) 


wtere  the  single  subscript  now  denotes  a  specific  point  in  the  structure. 


Similarly, 


*11  *12  1  T*  31  *  32  1  1  f*33  * 


-  -U  b- ■  JU41  *4J  U13  bU 


=  - 


*21  *22  *  31  *32  1  *33  *  34 
*31  *22  L^41  *42  **J  *44 


(Hll) 


(HI  2) 


where  *,]  are  3x3  subraatrices  cf  B'  =  X;'’2X2X1X2 

To  obtain  the  frequency  response  functions  far  stations  4,  5  and  6,  transforma ticn 
matrices  are  utilized.  First  consider  the  farces  at  junction  1  as  shown  in  Fig.  H2.  Frcm  equilibrium 
considerations, 


(HI  3) 


Similarly,  for  junction  1’, 


..i  ------  v  .  -•  -■  -•  , 


wtere  tie  subscripts  R  and  L  stand  for  right  and  left,  respectively,  cf  pants  1  and  1*.  Eqns.  (HI 3) 
and  (H14)  can  be  combined  as 


ipi)  _'Pi\  .Pi\ 

Ipi  Jl  lPl,  K  {Pi,  L 


(HI  5) 


The  state  vectors  to  the  left  and  right  of  junctions  1  and  1*  are  related  by  the  transfer  matrix  X2 
such  that 


\ 

d[ 

p[ 

Pi 

* 

/  0  0  o]  di 

0  I  0  0  ^ 

GiC{xGi  -GiC{lC  iG3  1  0  p\ 


(HI  6) 


from  which 


' 

Pi 

Pi 

v  > 

X 

GiC{lG\  -Gt,C2xC  \Gi  /0  d\ 

G£tCilGi  C«(Cj-C«Cr1C1)C30  I  p[ 


(H 17) 


Combining  eqns.  (Hi 5)  and  (HIT), 


G  ,C  {XG  i  -GtCiH}}  0  0  dl 

C:CtCfxG l  G4(C3-C£{'Ci)G30  0  p[ 


Since  displacements  are  continuous  in  crossing  over  junctions  1  and  i ’, 
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Eqns.  (HI 8)  and  (HI 9)  can  be  combined  as 


di 

/  0  0  o' 

di 

dv 

Pi- 

'  = 

0  /  0  0 

G<CilGl  -GtC^Cfi  J  0  0 

dy 

Pi- 

Pi 

> 

G£<C{lGi  Gi(C3-C4Ci-1C1)Gj°  0 

Pi 

Keeping  in  mind  that  1  corresponds  to  member  I  and  1*  corresponds  to  member  n,  eqn. 
be  rewritten  as 


d, 

I  0  0  0 

f  \ 

d, 

d„ 

0  I  0  0 

da 

Pu 

G<C{lGl  -GtC{1ClGi  0  0 

Pu 

Pi 

J 

G££i'G,  Gi(C j-C4C ilC i)Gj  0  0 

Pi 

<  )  L 

Now  let  the  transfer  matrix  in  eqn.  (H21)  be  X2  for  station  1-1', 


r  ' 

d, 

\ 

d, 

dtt 

=  *2*1*2 

dn 

Pu 

Pu 

Pi 

> 

i-r 

Pi 

(HI  9) 


(H20) 

(H20)  can 


(H21) 


(H22) 


Let  b,)'  be  3x3  submatrices  c i  Bm  =  X 2X1X2,  then,  since  there  are  no  force  constraints  at  station  0- 

0’, 


38- 


Tte  directional  relationships  between  the  local  state  vectors  in  member  11*  to  that  ci  the 
global  ores  in  member  II  at  junction  point  1’  is  shown  in  Fig.  H3.  Inspection  d  Fig.  H3  gives  the 
following  relations, 

-vvr  =  -uv  Mv  =  -wr  ,<k-  -  *|V 

cr 


’  ^ 

u 

0  1  0 

u 

-Av 

■  = 

-10  0 

— tl) 

VT 

r 

0  0  1. 

<1 

\ 

and  Vv  =NV  ;M t.  =  -My  ;Ny  -  ~Vy 


r 

> 

M 

-1  0  0 

M 

V 

•  = 

o 

o 

V 

s 

r 

[o  -1  oj 

.V 

Combining  eqns.  (H24)  and  (H25) , 


A 

O 

o 

'd 

P)v 

[o  c;j 

P, 

(H24) 


(H25) 


CH25) 


where 


C[  = 


0  1  0 

-10  0 
0  0  1 


and  where 


r.-  v  '.*v  •.*  v 


y~.  • '  ■"  t’ i  '  ».■  »■’ «.■  i  ■. '  ^'w.1 1»  ti«j  i  ,  ■■  ..t  .....  .j  .y 
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Gi  = 


-10  0 
0  0  1 
0-10 


Referring  back  to  Fig.  H3,  at  station  5,  which  involves  the  transfer  cf  state  vectors  in 
station  1’  through  a  length  cf  €/2  fran  1’  to  5, 


0,  -rwf3.. 


or 


£).-«  Cl 


(H27) 


wtere  Tm  is  the  transfer  matrix  for  a  Timoshenko  beam  cf  length  €/2.  Since  frcm  eqn.  (H23) , 


• 

j  ^21*22  U 

U 

\p)v  = 

Substitution  af  eqn.  (H28)  into  eqn.  (H27)  gives 


(H2S) 


Gi  o 

*31*32 

-1 

• 

^33^34 

fill 

.  r  =  -rw 
v>Js 

o  g2 

*4i*42 

p43^«4 

c~c 


fran  which 


=  -rw 


c; 

0 

*31*32 

-i 

b-abu 

0 

Gj 

f)3>H 

k  4i  *  *2 

p43^44 

(H29) 


(H30) 


Similarly,  for  station  4, 


-140- 


=  -r1-7 


Gj  o 
O'  Gi 


MT  bZ 

h r 


*  31  *41 

*32  *42 


-1 


&33  M 

^44  ^44 


(H31) 


where  btf  are  3x3  submatrices  cf  By  ~  X2.  Far  station  6, 


h%  ^'(w)  =  -r1*7 


o 

o 

_ i 

Ml  Ml 

*31  *32 

-1 

633  fcj4 

[0  g2. 

.Mi  M*. 

Ml  *42 

&43  bu 

where  b *  are  3x3  submatrices  c i  Bv  =  XiXiX2Xx. 


(H32) 


To  obtain  the  state  vector  at  locations  to  the  right  cf  excitation  station.  C-C\  it  is  more 
convenient  to  use  transfer  matrices  from  right  to  left.  Since  the  derivation  follows  the  same 
procedure  outlined  previously,  it  is  not  discussed  further  here. 
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Fig.  H 2  Forces  at  location  1. 
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APPENDIX  1 

TRANSFER  MATRIX  FOR  3-D  WAVE 
PROPAGATION  IN  A  TETRAHEDRAL  TRUSS 

The  transfer  matrices  for  wave  propagation  in  a  tetrahedral  truss  are  derived  in  this 
appendix.  Fig.  II  shows  a  tetrahedral  truss  with  three  repeating  periodic  units,  which  is  to  be  used 
in  the  derivation  of  the  transfer  matrices.  The  tetrahedral  truss  is  assumed  to  be  made  up  of  identi¬ 
cal  elastic  bars  with  distributed  mass  and  circular  a  css  sections,  each  cf  length  i.  Each  connecting 
bar  has  modulus  of  elasticity  E ,  mass  density  p,  shear  modulus  G,  cross-sectional  area  A,  second 
moment  of  area  inertia  about  the  Xor  y  axis  J ,  second  moment  c£  area  inertia  about  the  Xaxis  JT , 
and  radius  of  gyration  about  the  x  or  z  axis  i .  The  crass-sectional  dimension  of  each  connecting  bar 
is  assumed  to  be  small  compared  with  its  length.  Fig.  12  shows  the  sign  conventions  for  the  forces 
and  displacements  in  a  connecting  bar  in  the  global  xyz  coordinates.  Each  bar  has  longitudinal  dis¬ 
placement  u,  transverse  displacements  v  and  w  in  the  y  and  z  directions,  respectively,  and  angles  of 
rotation.  4,  4  and  0  about  the  x,  y  and  z  axis ,  respectively.  In  addition,  each  bar  also  has  moments 
My  and  M,  about  the  y  and  z  axes,  respectively,  torsion  T ,  axial  force  N  and  shear  farces  Vy  and 
V ,  along  the  v  and  z  axes,  respectively. 

In  the  analysis  of  wave  propagation  in  a  tetrahedral  truss  (Fig.  II),  three  transfer 
matrices,  namely,  V ; ,  V :  and  V3  are  developed.  Fig.  13  shows  an  arbitrary  segment  cf  the 
tetrahedral  truss  and  Fig.  14  shows  the  segment  which  is  sectioned  along  planes  parallel  to  the  yz 
planes  into  the  substructures  responsible  for  the  transfer  matrices  Vu  V2  and  V2.  Referring  to  Fig. 
14,  the  sections  are  made  by  cutting  the  periodic  unit  (Fig.  13)  along  planes  GAnlDx ,  ERKRFRLR 
(or  ElKlFlLl)  and  HBLJCL  which  are  all  parallel  to  the  yz  planes.  The  subscripts  L  and  R  used 
for  points  E ,  A,  D.F,X  ,L,B  and  C  denote  points  just  left  and  right  of  these  points,  respectively. 
Thus,  the  first  transfer  matrix  V;  represens  the  txxnsfer  cf  state  vectors  in  members  which  join 


■:  v  1/  J  I  1 .1  L>»  !■ 


€ 
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members  I  and  HI.  For  example,  members  AD  and  BC  each  result  in  a  transfer  matrix  Vt.  The 
second  transfer  matrix  V2  involves  the  transfer  d  state  vectors  in  four  bars,  each  of  length  i/2  in 
members  I  through  IV  together  with  members  which  connect  member  I  to  member  II,  member  HI 
to  member  H,  member  I  to  member  IV  and  member  EH  to  member  IV.  Referring  to  Fig.  14,  this 
includes  members  AkKl,  1Fl,  DmLl,  GiEl,  A„Fl,  DkFl,  AxEc,  DkEl.  Lastly,  transfer  matrix 
V3  is  responsible  for  the  transfer  of  state  vectors  in  four  bars,  each  cf  length  i/2  in  members  I  and 
HI  together  with  members  which  connect  member  H  to  member  I,  member  El  to  number  EH, 
member  IV  to  member  I  and  member  IV  to  member  EH.  For  example,  referring  again  to  Fig.  14, 
this  includes  members  KkBl,F„J  ,LkCl,Ei,H  ,FkBl,FkCl,  EkBl  and  EKCL . 


f  .V 


transfer  matrix  v, 


For  the  flexural  and  longitudinal  3-D  vibration  in  an  elastic  bar,  corresponding  to  a  state 


vector  z  where  z 


and 


ft/1 
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M, 
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.  P  = 

V, 

6 

9 

U  j 

the  transfer  matrix  T  is  given  as  [5]. 
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1  _  1^.  »  .  B  .  « 
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T  = 


C, 

C3 


C2 

c4 


where  Cu  C2,  C3  and  C<  are  6x6  sub  matrices  such  that 
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a 
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0 
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0 

0 
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co-rcj  0  0  0  €[ci— (o*+-t)c3] 
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a  = 
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1  =  \/\/34+^-(ov-t)1  +  \ 
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X^+Xf 
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(o+r) 


Consider  member  AD  (Fig.  13)  which  is  representative  cf  transfer  matrix  V x.  Fig.  15 
shcr*s  the  directional  relationships  between  the  global  state  vectors  and  the  local  state  vectors. 


Inspection  cf  Fig.  15  gives  the  following  relations, 


u  —  -v  ;  v  =  u  i  -w  =  -w 
4>  =  -<j»  ;  =  <{>  -.  0  =  0 

M,  =  M,  ;My  =T  ;T  =  -M, 
V,  =V,  ;V,  =N  ;S  =  -V, 


In  matrix  form,  for  junction  D , 
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V, 
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Vy 

u 
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D 

Similarly,  for  junction  A, 


sphere  G;  -  G}  such  that 


G i  -  C3 


0  -1  0  0  0  0 

1  0  0  0  0  0 

0  0  1  0  0  0 

0  0  0  0  -1  0 

0  0  0  1  0  0 

.0  0  0  0  0  0 


and 


G  2  —  -G4 


1  0  0  0  0  0 

0  0  1  0  0  0 

0  -1  0  0  0  0 

0  0  0  1  0  0 

0  0  0  0  0  1 

0  0  0  0  -1  0 


For  member  AD , 


Ci 

Cl 


C2 

c4 


where  Ci,  Ci,  Ci  and  C4  are  defined  previously. 
From  eqn.  (112), 


dr)  —C\d,\  +  C  lP  4 
PO  ~  C  id  A  +  C  t,pA 


Multiplying  eqn.  (113)  by  C2_1 , 
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Pa  =C?dD  ~C?C,dA 


Substituting  eqn.  (115)  into  eqn.  (114), 


po  =C<C{ld0  +  (Cj-C.CftO^ 


Using  eqns.  (115),  (16)  and  eqns.  (18)  and  (17), 


Pa  —  Gt,C{xG\do  —  Gt,C{xCiG^iA 


Using  eqns.  (116)  and  (16)  and  eqns.  (18)  and  (19), 


Pd  ~GjCtCz^Gid0  +C2(C3  —  C«C2  xCi)Gt4a 


Eqns.  (117)  and  (118)  give  the  force  response  at  A  and  at  D  due  to  the  displacements  at  A  and  D. 

New  consider  Fig.  16,  which  shows  the  forces  at  junctions  A  and  D.  From  equilibrium 
considerations, 


Pa  1  (Pa  )  (Pa  ) 

Pa)  ^  [Po  )  ( Po)  1 


where  the  subscripts  L  and  R  denote  points  just  left  and  just  right  of  A  ct  D. 


Keeping  in  mind  that  A  is  in  member  I  and  D  is  in  member  m,  using  eqns.  (117),  (118) 


and  (119), 


-V 


*  * 
/V,. 
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0  0  0  I  0  0  Ol \p„ 

G2(C 3-C4C 2-1C j)C3  0  G^Cf'Gj  0/00  PJ,J 
=  n  a  a  00/0'- 


u  u  u  u  yj  1  u 

-GtCi'CxGi  0  G4Cj_1Gi  0  0  0  / 


where  the  subscripts  I,  H,  HI  and  IV  denote  member  numbers. 


In  crossing  the  joints  A  and  D,  since  displacements  are  continuous, 


/  0  0  0 
0/00 
0  0/0 
0  0  0  / 


Combining  eqns.  (120)  and  (121), 


Pm  j  R 


G:(C3-C<C2-1C1)G3 

0 

-GiC{'CiG) 

G2(C3-C4C2-lC1)G3 


0  0  0  0  0  0  j  d, 

0  0  0  0  0  0  '  d„ 

0  0  0  0  0  0;  d,„ 

I  1  0  0  0  0  '  dy 

0  0  I  0  0  0  i  py 

^1  0  0  1  0  0  |  pm 

0  0  0  0  I  0  ,  pa 

GjCf'Gi  0  0  0  0  1  j  p, 

0  0  0  I  0  0  !  p„, 


Thus ,  from  eqn.  (122), 
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I  0 

0  I 

0  0 

Vi  =  0  0 

0  0 

G2(G3-C«C2-1C1)G3  0 

0  0 

-GtCi'CiGi  0 


0 

0 

I 

0 

0 

G2C4C2-1G1 

0 

GtC{lGi 


0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
1  0  0  0  0 
0  10  0  0 
0  0  10  0 
0  0  0  1  0 
0  0  0  0  1 


TRANSFER  MATRIX  V2 


(123) 


Fig.  17  show  a  substructure  which  is  representative  cf  transfer  matrix  V2.  Consider  first 
memebers  AkKl  ,  IFL ,  DmLl  and  GEL .  Ncxice  that  these  members  are  parts  cf  main  members  I, 
II,  in  and  IV;  respectively,  each  cf  length  111. 

From  eqn.  (II),  the  transfer  matrix  T  for  longitudinal  and  flexural  vibration  in  a  bar  at 
lenth  /  is  given  by 


Ci  C2 
C3  C« 


where  the  submatrices  Cj,  C2,  C3  and  C«  where  defined  previously.  Thus,  for  a  bar  cf  length  tl 2, 
the  transfer  matrix  T '  is  given  by 


Cj  Cj 

c3  Ci 


(124) 


where  the  submatrices  Ci,Cj,C3  and  C«  are  obtained  by  substituting  i/2  far  (  in  submatrioes 
C;>  Ci,  C3  and  C4,  respectively. 
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Using  eqn.  (124),  the  foUowing  relationships  are  obtained, 

ft.- 

ft- 

d) 

phi 

Keeping  in  mind  that  AM  ,1 ,  D*  and  G  are  all  left-end  points  of  the  substructure  (shown 
in  Fig.  17)  in  numbers  I,  n,  HI  and  IV,  respectively  and  KL  ,  FL  ,  Ll  and  EL  are  all  right  end  points 
cf  the  substructure  in  members  I,  II,  HI  and  IV,  respectively,  and  using  eqns.  (125)  to  (128). 

Cj  0  0  0  0  0  0  Ci*1  f  d, 

0  C|  0  0  0  0  Ci  0  d„ 

0  0  C[  0  0  Cj  0  0  d„, 

0  0  0  ci  C:  0  0  0  </yv 

0  0  0  d  ci  0  0  0  Spv 

0  0  Cj  0  0  C't  0  0  p,„ 

0  Cj  0  0  0  0  C«0|p// 

cj  o  o  o  o  o  o  c;i  P,  I 

v  J 

Ml  (129) 


where  the  subscripts  Mm  and  ML  denote  the  right-end  and  left-end  points  of  the  main  members 
aniy. 


Now  consider  members  /t#  FL  ,  FL  ,  E  and  DK  EL  ,  which  are  the  remaining 


members  to  be  analyzed.  Notice  that  these  members,  each  cf  length  l ,  join  the  main  members 
together  and  the  local  state  vector  coordinate  cf  each  does  net  coincide  with  the  global  state  vector 
coordinates  in  the  four  main  members.  Thus,  with  respect  to  each  number,  transformation  matrices 
have  to  be  utilized  in  order  to  consider  the  transfer  a i  state  vectors  from  one  paint  in  a  main 
member  to  another  point  in  another  main  member  (fer  example,  from  As  in  main  member  I  to  FL 
in  main  member  II  through  member  AM  FL ). 

Consider  first  member  AR  FL  .  Fig.  18  shows  the  orientation  of  member  AMFL  in  the  glo¬ 
bal  xyz  coordinate.  The  directional  relationships  between  the  state  vectors  in  the  local  coordinates 
xyz  and  in  the  global  coordinates  xyz  at  fL  are  shown  in  Fig.  19.  Notice  that  the  local  xyz  axes 
are  chosen  such  that  the  y  axis  is  in  the  yz  plane.  To  find  the  transformation  matrices  needed  to 
convert  the  state  vectors  from  the  local  coordinates  to  these  cf  the  global  ones,  first  assume  that  the 
xyz  coordinates  is  rotated  about  the  x-axis  to  a  new  x'y'z'  coordinates  such  that  now  the  new  y  -axis 
coincides  with  the  local  y  axis.  Notice  that  member  AK  FL  is  now  in  the  x  y ' plane.  Fig.  110  shows 
the  directional  relations  between  the  local  and  the  rotated  global  state  vectors  at  FL  .  Examination 
of  Fig.  110  gives  the  following  relations: 
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u  =  (cos  60°)  (u)  +  (-ccs  30°)(-w‘) 
a  u  =  (y)0O  +  (~)(-w)  ,  (130) 

~-W  =  (cos  30°) («")  +  (-ccs  60°)(-w‘) 

.  j— 

or  -W  =  (~)(u)  +  (|)(^v)  ,  (131) 

4>  =  (cos  60°)  (4>)  +  (-ccs  30°)  (-S') 

or  4>  =  ( j)(40  +  (-y-)(-«)  ,  (132) 

0  =  (-ccs  30°)  (0)  +  (cos  600)  (-S') 

°r  8  =  (~)(0  +  (y)(-9)  , 


(133) 


and  Mi  =  (cos  60 °)(A/,)  +  (-ccs  30°)  (7) 


or  Mi  =  (j)(A/,)  +  (~)(T)  , 

7  =  (-ccs  2CP)(Mt)  +  (cos  60°)  (7) 
or  r’  =  (-~-)(A/,)  +  (1)(7)  , 
v;  =  (cos  60°)(v\)  +  (-<cs  30°)(7) 
cr  v;  =  (|)(V,)  +  (- y-)(T )  , 

N'  =  (-ccs30°)(V,)  +  (cos  60°)(7) 
cr  AT  -  (^-)(Vf)  +  (y)(7)  , 

Using  eqns.  (130)  to  (133), 


i 

0 

^3 

0 

0 

0 

/  \ 

2 

2 

• 

u 

0 

1 

0 

0 

0 

0 

V 

0 

1_ 

0 

0 

0 

— w 

2 

2 

6 

0 

0 

0 

_1_ 

0 

^3 

2 

0 

0 

0 

0 

1 

0 

k  0  . 

n 

0 

0 

0 

V3 

0 

j_ 

. 

2 

2 

and  using  eqns.  (134)  to  (137), 


1  o  -22. 

2  2 

0  1  0 

4o  i 

2  2 

0  0  0 

0  0  0 

0  0  0 


0  0  0 
0  0  0 
0  0  0 


0  T 


To  obtain  tie  relationships  between  the  unrelated  global  state  vectors  and  the  local  state  vectors,  the 
x'y'z'  coordinates  are  rotated  back  to  the  original  xyz  position.  This  can  be  done  using  another  set 
cf  transformation  matrices.  Fig.  Ill  shows  the  relationship  between  the  global  state  vectors  and  the 
rotated  state  vectors.  As  before,  examination  of  Fig.  110  gives 

1  0  0  0  0  0 
t  V  o  V?  0  0  0  ,  V 


Y3 


8  )rL 


0  0 


0 

0 

0 

Vf 0 

0 

0 

0  1 

0 

0 

0  0 

1 

v3 

Vf 

0  0 

-Vf 

1_ 

v3 
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and  where 


2  v  r  tr  \2 

v3  „  1 


0  0 


G  2 - Gt  — 


M  i 

■iVj  Vf  -f 


0  f 


From  eqn.  (II),  in  the  local  xyz  coordinate, 


d)  _  [C,  C2](d) 

p)rL  [C3  Ci\  \p)A 


Thus,  from  eqn.  (150), 


dfL  =  CkdA„  +  CzpA„ 
and  prL  =  CydAn  +CipAlt 


Multiplying  eqn.  (151)  by  C2-1 , 


P.kh  ~  C  2  lJfL  ~C2lC  i  dAft 


Substituting  eqn.  (153)  into  eqn.  (152), 


*• 


Pfl  =C<C{xdFL  +  (C,  ~  C *C i'C i)d A„ 

(154) 

Using  eqns.  (153),  (146),  (148)  and  (149) 

Pak  =GAdxG\drL  —  G4CiCiG3(fA* 

(155) 

Using  eqns.  (154),  (146),  (148)  and  (145), 

PrL  —  G^C ACixG\dfL  +  G2(C3  -C4Cf1C1)G3d/tJ, 

(156) 

If  the  same  procedure  is  carried  out  fcr  members  DKFL,  AREL  and  DKEL,  and  if  far 

tie  analysis  of  each  member,  taking  the  left  end  point  and  the  origin  cf  the  local  coordinate  and  the 

direction  along  the  length  c£  the  member  as  the  local  i  direction,  with  the  y  axis  in  the  global  yz 

plane,  the  following  relations  are  obtained. 

For  member  DRFL , 

Po„  =G;CixG[dri  -G:CixCfildDK 

Pfl  -  GiC aC ixG[dFL  +  G (C 3  -  C.C^C-JG^ 

057) 

(158) 

T 


$ 


T 


cr  =g;= 


0  Vf  i  0  0 

T  I V?  "iVy  o  o 

n  n  ft  11 


0 


2 


and  where 


-iVf  v?  i 

-M  vf  i 

_^5  n  1 


0  0  0 


t 


For  member  £>*  , 


0  0  0 
0  0  0 


0  0 

-IV?  -fe 

-IVIVf 


Pon  =Gl'CixG['dtL  -  Gl'Ci^C  \G”dDn 

Pn  —  G  ”C  iC  2~lG  l"d  El  +  G2"CC3  —  C4C  2_1C  i)G"d0/l 


where 
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and  where  G"  =  ~0" 


Eqrs.  (155)  to  (162)  give  the  relations  between  the  induced  forces  at  the  end  points  cf 
member  ARFL,DRFL,  AnEL  and  DREL  due  to  displacements  at  AR  and  Fl,Df  and  FL,  AR  and 
El,Dr  and  EL,  respectively.  Keeping  in  mind  that  AK  refers  to  a  point  in  main  member  I;  DR 
refers  to  a  point  in  main  member  HI,  EL  refers  to  a  point  in  main  member  IV,  and  FL  refers  to  a 
point  in  main  member  II,  and  introducing  a  set  cf  matrices  D\  through  D16  such  that 


V 


G 
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;c2~ig; 


■■  -<JiCilC:Gi 


GjC^Cj-'g; 


G}(Ci  —CtC{1Ci)Gi 
■ GlCfG i 
■  -g;c2-ic1g; 


■g;c<c2-ig; 

>g;(Cj  -C4Ci-lCl)G,* 


=  g4"c2-£g 


w 

1 


=  -g;cj-ic1g3- 


^G'XtC^G" 


=  G2'(C3  —  C4C2_1C  i)G3’ 
=c;'c2-1Gf' 


=  -GrWxG,- 


=  G2"C4C2-’Gr 


(C3  -  C4C2_1Ci)G3",  eqns.  (155)  through  (164)  can  be  expressed  in  matrix  farm  as 


0 

0 

0 

0 

0 

0 

0 

0 

d, 

0 

0 

0 

0 

0 

0 

0 

0 

d„ 

0 

0 

0 

0 

0 

0 

0 

0 

dm 

0 

0 

0 

0 

0 

0 

0 

0 

\d„ 

’  _ 

Da 

0 

016 

011 +01J 

0 

0 

0 

0 

Prv 

0 

0s 

06+014 

013 

0 

0 

0 

0 

Pm 

Dt 

03+07 

0. 

0 

0 

0 

0 

0 

Pii 

Di 

0 

0, 

0 

0 

0 

0 

/ 

M 

where  the  subscripts  CL  and  CM  denote  the  left  and  right  ends  c f  the  connecting 
members  which  join  the  main  members  in  a  section  c£  a  periodic  unit  which  results  in  transfer 
matrix  V2.  New,  adding  the  contributions  frem  members  AkKl,IFl,DkLl  and  G:EL  using  eqn. 
(129),  to  that  in  eqn.  (163),  transfer  matrix  VL  can  be  obtained.  Notice  that  the  addition  cf  the  ele¬ 
ments  in  the  matrices  in  eqns.  (163)  and  (129)  is  in  effect  adding  the  contributions  from  the  consti¬ 
tutive  members  which  make  up  transfer  matrix  V  t.  Thus, 


d, ) 

* 

c; 

0 

0 

0 

0 

0 

0 

c2 : 

d„ 

0 

c; 

0 

0 

0 

OCj 

C: 

o  i 

dm 

0 

0 

c; 

0 

0 

C2 

0 

o  ! 

dp/ 

0 

0 

0 

c; 

c2 

0 

0 

0 

Prv 

012 

0 

016 

011+013+^3 

c« 

0 

0 

0 

Pm 

0 

03 

06+014 +Cj 

013 

0 

c; 

0 

0  :  ! 

Pn 

04 

3 

03 

0 

0 

0 

c4 

o  1 

P  ’ 

Di+D  io^C  3 

Di 

0 

02 

0 

0 

0 

c;.  . 

iron  which 


FJ r’VZ'T’y 
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V,  = 


c; 

0 

0 

0 

0 

0 

0 

Ci 

0 

cl 

0 

0 

0 

0 

Cl 

0 

0 

0 

Cl 

0 

0 

ci 

0 

0 

0 

0 

0 

Cl' 

Ci 

0 

0 

0 

£>12 

0 

016 

011+015 

■+C3Ct 

0 

0 

0 

0 

0J 

06+014^3 

0 13 

0 

c; 

0 

0 

04 

D3+D7+C3 

08 

0 

0 

0 

C4 

0 

02+010 

+Ca  Di 

0 

09 

0 

0 

0 

C4 

(164) 

where  the  subscripts  /?  and  L  denote  the  right  and  left  ends  cf  the  members  in  a  section  cf  a 
periodic  structure  which  results  in  transfer  matrix  V2. 

TRANSFER  MATRIX  V3 

Fig.  112  shows  a  substructure  which  is  representative  cf  the  transfer  matrix  which  is 
representative  of  transfer  matrix  V3.  The  derivation  cf  the  transfer  matrix  V3  is  similar  to  that  of 
the  transfer  matrix  V2.  Referring  to  Fig.  112,  members  KKBL,FgJ ,L„Cl  and  EKH,  which  are 
parts  cf  main  members  I,  II,  HI  and  IV,  respectively,  are  considered  first.  Members 
FrBl,F„Cl,EkBl  and  EKCL  are  then  considered.  The  transfer  matrix  V3  is  then  obtained  by 
adding  the  contributions  cf  all  eight  members,  as  was  done  for  transfer  matrix  V2.  Since  the  pro¬ 
cedure  is  similar  to  that  cf  the  previous  section,  it  will  not  be  repeated  here.  Transfer  matrix  V3  is 
obtained  as 


V3  = 


c; 

0 

0 

0 

£9 

0 

Ei 

E^E  K- 


0 

c; 

0 

0 

0 

£s 

Ei^Et-rCi 

.C 


0 

0 

cl 

0 

£u 

Er^E^-^C  3 


r, 

—  5 


0 

0 

0 

ci 


0  0  0  C: 


0 


0  0  C2  0 

0  c2  0  0 
c2  0  0  0 

EiQ-*-En-fC 3C 4  0  0  0 

£;s  0  C'<  0  0 

0  0  0  c\  0 

El  0  0  0  C 4 


(165) 
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£ .  =H<C{lH[ 

E2  =  -H'£ixCxHi 
£3  =  H£AC?Hl 
£4  =W2(C3  -C4C:-lC,)//i 


£5  =  tf  ;c2-iffi* 


£6  =  -h,v2-icxh; 
En  =h:c<c£1h; 


E»  =H: ;(Cj  -C4C2-1C1)//3" 


f—  H  -lu  "> 

9  —  «4  u  2  n  1 


£w  =  -hI’c^c.h; 


» 


£„  -h;c&?h: 

E-1  -ffj"(C3  -  CxC2xCx)H2 


£:J  =Hrc{'Hr 

£14  =  -H"C2-lClH" 
t 

EIS  =  H2C<C2-'H{- 

£-.6  -»r(C,  ~C,C£'-CX)H" 

> 


and  HI  through  Hi’  are  transformation  matrices  such  that 
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re  defined  in  the  previous  section  on  the  transfer  matrix  V 2. 
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APPENDIX  J 

UST  OF  COMPUTER  PROGRAM 

This  section  gives  a  Basic  language  computer  program,  named  PROG1,  BAS,  which 
calculate  the  frequency  response  functions  for  an  elastic  rod.  Fig.  J1  shows  a  schematic  cf  the  rod 
and  the  symbols  far  various  locations  used.  PROG1BAS  calculates  the  frequency  response 
functions,  according  to  discrete  frequency  steps,  far  locations  A,B ,C ,  D,  EL,  EK  ,F  and  G,  where 
El  and  Er  are  points  just  left  cf  E  and  just  right  of  E,  respectively.  The  values  of  the  frequency 
response  functions  are  stored,  according  to  frequency,  onto  files,  one  file  for  each  location.  The 
letter  before  the  decimal  point  in  each  file  name  denotes  the  response  location  and  the  ward  used 
after  the  decimal  paint  denotes  the  type  af  response  which  is  being  stored.  If  a  file  name  ends  with 
the  word  one,  the  response  cf  the  longitudinal  displacement  is  stored  whereas,  if  a  fiile  name  ends 
with  the  word  two,  the  response  cf  the  axial  force  is  stared. 


PROG1.BAS 


on  error  goto  84 
defdbl  t , e , w , m , a , b , 1 
def int  i  ,  j 
option  base  1 
din  t(2 ,2) 
read  e  ,o  ,  1 ,  a 

data  10.8e6,2.6e-4,10,9.75e-2 
input  "Enter  the  starting  point 
input  "Errter  the  terminating  point 
10  for  k2=k1  to  k3  step  10 
w=k2*1 00 
for  i=1  to  2 
for  j  =  1  to  2 
t( i *  j  )=0 
next  j 
next  i 
1=1*2 

gosub  1000 
bl  =  t(2 ,2) 

1= 1/2*6 
gosub  1000 
a1=t(2 ,2) 

open  "a",  1,  "freqa.one" 
write  #1 ,  w ,  -bl  / a  1 
close  1 
1=1/6 

gosub  1000 

open  "a",  1,  "freqb.one" 
write  H 1,  w,-bl/a1*t(l ,2) 
close  1 

open  "a",  1,  "freqb.two" 
write  #1,  w,-bl/a1*t(2,2) 
close  1 
1  =  1*2 

gosub  1000 

open  "a",  1,  "freqc.one" 
write  #1 ,  w , -bl /a  1  * t( 1  ,  2 ) 
close  1 

open  "a",  1,  "freqc.two" 
write  if  1,  w,-bl/a1*t(2,2) 
close  1 
1= 1/2*3 
gosub  1000 

open  "a”,  1,  "freqd.one" 
write  ft  1  ,  w,-bl/a1*t(l,2) 
close  1 

open  "a",  1,  "freqd.two" 

write  f!  1,  w ,  -  bl  /  a  1  *  t  (  2 , 2  ) 
close  1 
1=  1/3*4 

POSUb  1000 


open  "a",  1,  "freqel.one" 
write  #1  ,  w,-bl/a1*t(l,2) 
close  1 

open  "a",  1,  "freqel.two" 
write  #1,  w , -bl/al * t( 2 , 2) 
close  1 
b2=- t( 1,2) 

1= 1/4*6 
gosub  1000 
'  a2  =  t( 1 , 1  ) 


open 

"  a  "  , 

1 ,  " f req  g . one" 

write 

#1  , 

w, b2/a2 

close 

1 

1=1/6 

gosub 

1000 

open 

II  a  II 
a  > 

1,  "freqf.one" 

write 

w,b2/a2*t(1 ,1 ) 

close 

1 

open 

It  a  *• 

t 

1 ,  " f reqf . two" 

write 

w,-b2/a2*t(2 ,1 

close 

1 

1  =  1*2 

gosub 

1000 

open 

it  a  tt 
®  » 

1 ,  "f reqer . one 

write 

#1  . 

w,b2/a2*t(1 ,1 ) 

close 

1 

open 

■;a  ”  , 

1 ,  " f reqer . two 

w  ri  te 

w,-b2/a2*t(2, 1 

close 

1 

1  =  1/2 

next 

k2 

goto 

86 

84  kl =k2+  1  0 

resume  10 
86  end 

1000  th=l*w*sqr(m*a/e) 
t( 1 , 1 )=cos( th) 
t(1  ,2)=l*sin(th)/ e/a/th 
t(2,1  )=-m*l*w*w*sin(th)/'  t  J 
t  (  2 , 2  )  =  t  (  1  , 1  ) 
return 
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APPENDIX  K 

SOME  PROPERTIES  OF  TRANSFER  MATRICES 

CROSS- SYMMETRY 

If  a  segment  of  an  element  is  symmetric  about  a  plane  at  its  midlength  as  shown 
schematically  in  Fig.  Kl,  then  it  is  always  possible  to  obtain  a  cress -symmetric  transfer  matrix, 
which  represents  the  element,  by  a  suitable  ordering  and  sign  convention  of  the  components  of  the 
state  vector  [7].  The  transfer  matrix  is  cross-symmetric  in  that  its  elements  are  symmetric  about  its 
crass -diagonal.  Thus,  if  the  elements  in  a  crass -symmetric  transfer  matrix  T  are  r(J,  where 
i  =1,2,.../j  for  the  rows  and  j=l2,—P  for  the  columns  and  n  is  the  dimension  of  the  transfer  matrix 
(transfer  matrices  are  square  matrices),  then 

lij  =  (XI) 

For  example,  the  transfer  matrix  for  an  elastic  rod  can  be  observed  to  be  cross- 
symmetric  using  the  state  vector  as  defined.  Similarly,  it  can  be  shown  that  the  transfer  matrix  for 
a  Timoshenko  beam  is  cross-symmetric  with  a  suitable  ordering  cf  the  components  ctf  the  state  vec¬ 
tor. 

INVERSION 


For  a  given  element,  there  exists  two  possible  transfer  matrices  to  relate  either  (1)  the 
state  vector  at  the  right  end  to  the  state  vector  at  the  left  end,  or  (2)  the  state  vector  at  the  left  end 
to  the  state  vector  at  the  right  end.  Let  T  be  a  transfer  matrix  which  relates  the  state  vector  at  the 
right  end  to  the  state  vector  at  the  left  end  d  an  element.  Thus, 
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=  Tzl  (K2) 

wtere  zM  anj  zL  denote  the  state  sector  at  the  right  end  and  the  state  vector  at  the  left  end  of  the 
element,  respectively.  Eqn.  (K2)  can  be  rearranged  to  yield 

zt  =r-1z,  (X 3) 

Eqn.  (K3)  shows  that  the  state  vector  at  the  left  end  is  related  to  the  state  vector  at  the  right  end  by 
the  inverse  of  the  transfer  matrix  T . 

For  example,  the  transfer  matrix  relating  the  state  vector  at  the  left  end  to  the  right  end 
of  an  elastic  rod  is  the  inverse  of  the  transfer  matrix  relating  the  state  vector  at  the  right  end  to  the 
left  end  cf  the  same  rod. 

It  may  also  be  shown  that  transfer  matrices  are  non-singular  and  that  the  inverse  cf  a 
transfer  matrix  is  equal  to  its  adjoint  [7], 

VALUE  OF  DETERMINANT 

It  can  be  shown  that  the  determinant  cf  any  transfer  matrix  is  equal  to  unity  [7]. 

For  example,  the  determinant  cf  the  transfer  matrix  for  the  longitudinal  vibration  cf  an 
elastic  rod  (or  the  inverse  cf  the  transfer  matrix)  can  be  shown  to  be  unity  algebraically. 

Because  digital  computations  are  necessary  in  applying  the  transfer  matrix  method,  the 
numerical  evaluation  cf  the  determinant  cf  the  transfer  matrix  provides  an  opportunity  to  determine 
the  accuracy  cf  computer-aided  numerical  results.  Table  K1  shows  the  listing  of  a  Basic  language 
computer  program  named  DET1BAS,  which  calculates  the  determinant  of  the  transfer  matrix  far 
the  longitudinal  vibration  cf  an  elastic  rod  at  various  frequencies.  Table  K2  shows  the  listing  cf  a 
Basic  language  computer  program  named  DET2.BAS ,  which  calculates  the  determinant  cf  the 
transfer  matrix  for  the  flexural  vibration  of  a  Timoshenko  beam  at  various  frequencies.  The 


'As' 
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programs  are  evaluated  using  double  precision  on  all  variables  and  calculations;  that  is,  all  variables 
are!  calculations  are  carried  out  to  16  significant  figures.  Table  K3  shows  the  computer  values  cf 
the  determinant  far  the  rod  and  the  Timoshenko  beam  at  various  frequencies.  From  Table  K3,  it  is 
observed  that  the  determinant  of  the  rod  is  always  calculated  correctly  as  unity.  However,  it  is 
observed  that  the  determinant  of  the  Timoshenko  beam  is  calculated  to  be  unity  only  at  low  fre¬ 
quencies  such  as  less  than  approximately  6000  rad/sec  (approximately  1000  Hz). 

eigenvalues 

It  can  be  shown  that  the  eigenvalues  (that  is,  the  natural  frequencies)  corresponding  to 
vibrations  cf  a  structure  evaluated  using  the  transfer  matrices  relating  (1)  the  state  vector  at  the 
right  end  to  the  state  vector  at  the  left  end,  and  (2)  the  state  vector  at  the  left  end  to  the  state  vec¬ 
tor  at  the  right  end,  axe  identical  [7]. 

For  example,  the  natural  frequencies  cf  longitudinal  vibration  in  an  elastic  rod  shown  in 
Fig.  K2  can  be  evaluated  using  the  transfer  matrices  for  the  rod.  As  shown  in  Fig.  K2,  the  rod  is 
clamped  at  tte  left  erri  and  is  free  at  the  right  end.  The  state  vector  cf  the  rod  at  the  right  end  is 
related  to  the  state  vector  at  the  left  end  by  a  transfer  matrix  as 


cos  9 


— pfa? 


sin  9 


l  sin  9 
EA  9 

ccs  9 


u 

l-VJt 


(X4) 


where  u  is  the  longitudinal  displacement  N  is  the  axial  farce  and  the  subscripts  R  and  L  denote  the 
right  end  and  the  left  end  cf  the  rod,  respectively.  Imposing  the  boundary  conditions  ctf 


substitution  ctf  eqn.  (K12)  into  eqn.  (KI1)  gives 
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*<A/J  =  l^1)?  for  n= 12  (K13) 

From  eqn.  (K13),  the  natural  frequencies  for  longitudinal  vibration  d  an  elastic  rod  is 


for  n  =1 


(K14) 


Similarly,  the  state  vector  of  the  rod  at  the  left  end  is  related  to  the  state  vector  at  the 
right  end  by  a  transfer  matrix  as 


cos  8 
B  i  sin  8 

p*a?— 


I  sin  8 
EA  8 

cos  8 


Imposing  the  boundary  conditions  in  eqns.  (K5)  and  (K6)  gives 


0 

\n 


ccs  9 

# < 


(  sin  8 
EA  8 

ccs  9 


(K15) 


(XI 6) 


Eqn.  (KI6)  can  be  rewritten  as  the  foil  wing  two  equations 


0  =  cos  euK 

...  ,  -.sin  8 

and  NL  =  yfu*— — -u. 


From  eqn.  (X17),  nontrivial  sduticms  require  that 


(K17) 

(K1S) 


cos  9  =  0 


(K19) 


Since  eqns.  (KIO)  and  (K19)  are  identical,  the  same  natural  frequencies  as  given  in  eqn.  (K14)  will 
result.  Thus,  it  is  demonstrated  that  the  same  eigenvalues  (or  natural  frequencies)  are  obtained 
using  the  transfer  matrices  relating  either  (1)  the  state  vector  at  the  right  end  to  the  state  vector  at 
the  left  end,  or  (2)  the  state  vector  at  the  left  end  to  the  state  vector  at  the  right  end. 


v  .  Sv.‘I'-\-w*VVW"v 

1_  •  -  ^  .  «  .  .  «  .  /.  .*  .»  .»  *  ' 


.  *vwv: 


TabJeKl.  Computer  listing  cf  pregmm  DET1  BAS 


d  e  f  i  n  t  i,j 

defdbl  a 

option  base  1 

din  a  <  2,2) 

road  t, ar.ro,  l 

data  10e?  .  1  .  B4e-1  ,  3e-4 , 9  . 75 

input  "The  frequency  of  vibration  is  "  ;v»  :  print 

lprint  "The  frequency  of  vibration  is  *'  ;w  iprint 

for  i  *  1  to  2 

for  j  *  1  to  2 

a  < i  .  j  >  »0 

next  j 

next  i 

go  sub  1000 

det-atl ,l)«a<2,2>-a(l , 2 ) *  a ( 2 ,  1  > 

lprint  "The  determinant  of  transfer  matrix  is  ";det 
end 

1  000  t  h=*  1  *w*  s  q  r  (  r  o  /  e  ) 
a(  1  ,  1  )*cost  th) 
a  <  2  ,  2 ) «a ( 1  , l > 
a(l,2)«We/ar/th*sin(th) 
a  < 2  ,  1 >  * -e  *  a  r  *  th/ 1  *  s i n ( th) 
return 


Iprint 
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Tahlc  K2.  Computer  listing  d  program  DETZBAS 
deflnt 

def dbl  a , b , c , e  ,  g  ,  r 
option  base  1 
dim  a( 12, 12) ,b(l2,12) 
def  fncsh( f )  =  ( ex p( f ) 4  ex p( -f ) ) /2 
def  f nsnh( h )« ( ex p( h) -ex p( -h ) ) /2 
read  e , g , ar , bi , ro , r , 1 , ga 

data  10e7,3.84e7,1 .84e-1 , 2 . 85 e-3 , 3e-4 , 1 .  24  e-1  ,  9 . 75 , 5 . 92  e6 
input  "Enter  the  frequency  of  vibration  :  ";w 
lprint  "The  frequency  of  vibration  is  »  ";w 
gosub  2000 
for  i=1  to  12 
for  j=1  to  12 
a( i  ,  j )«0 
b(i,  i  )«0 
next  j 
next  i 
gosub  3000 
for  i*1  to  4 
for  j«1  to  4 
b(i,o)*a(i+4,  j+4) 
next  j 
next  i 

lprint  "The  beam  transfer  oatrix  is  " 

for  i= 1  to  4 

lprint  b(i,l);b(it2);b(i,3);b(i,4) 
next  i 

d1«(b(1,l)*b(2,2)-b(1.2)*b(2,l))*(b(4,4)*b(3,3)-b(3,4)*b(4,3)) 
d2-(b(l  ,1  )*b(2,3)-b(2f1  )*b(l  ,  3  ) )  *(  b  (  3 ,2  )*  b  ( 4 , 4  ) -b  (  4 , 2  )*  b  (  3 , 4  ) ) 
d3«(b(1,1)*b(2,4)-b(2f1)*b(lf4))*(b(3.2)*b(4f3)-b(4,2)*b(3.3)) 
d4«(b(l  ,2)*b(2,3)-b(2.2)*b(l  ,3)  >*(b(3.1  )*b(  4  ,4  )-b  (  4  , 1  )*b(  3 ,4  ) ) 
d5«(b(1,2)*b(2,4)-b(2,2)*b(l,4))*(b(3,l)#b(4,3)-b(4,l)*b(3,3)) 
d6  =  ( b  ( 1  ,3)*b(2,4)-b(2,3)*b(l  , 4  )  )  *(  b  (  3 , 1  )*b(4,2)-b(4,1  )*b  (  3 , 2  ) ) 
detsd1-d2+d34d4-d5+d6 

lprint  "The  determinant  of  the  beam  transfer  matrix  is  *  ";det 
end 

200C  suaro/ar 
al-( 1*2  )/e/bi 
b4=au* ( w‘2 ) * ( l*4)/e/bi 
s»su*(v*2)*(l*2)/ga 
t«su*(r‘2)*( w *2)*( 1*2) /e/bi 
th*l*w*sqr( to/ e ) 

11«aqr(sqr(b4-*'(((s-t)*2)/4))-(»4t)/2) 

12«  aqr( sqr( b4* ( ( ( s- 1 ) *2 )/4 ) ) ♦ ( S4  t )/2 ) 

10- 1 / ( ( 11 *2)4 (  12*2)  ) 

c0«l0*( ( 12*2)*  fncsn( 11 ) - ' 11 *2 ) *cos(  12 )  ) 
c1»l0*((  12*2  )/11*fnsnh(  11  )♦(  11  *2  )/12*sir.  i,  12  )  ; 
c2*lC*(fncsh( 11 )-cos( 12) ) 
c3=l0*(l/11*fnsnh(  11 )-1/12*sin( 12  )  ) 
return 
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Table  K2.  Cant.  Computer  listing  cf  program  DET2.BAS 

rem  subroutine  to  form  beam  transfer  matrix 
3000  for  i«1  to  12 
for  j*1  to  12 
a(i,j)-0 
next  j 
next  i 

a ( 1 , 1 ) scos ( th) 

•(4,4 )*«(1 ,1 ) 
a(9,9)*a(l ,1) 
a(l2,12)-a(1 ,1  ) 
a(1 ,12)»l*sin(th)/(e*ar*th) 
a(4,9)-a(1 ,12) 
a  ( 2 , 2 ) *  cO- 8* c2 

•  (  5 , 5 )«a ( 2 , 2 ) 
a(2,3)«l*(c1-(s+t)*c3) 
a ( 5  ,6  )«a ( 2 ,3  ) 

a ( 2 , 10 )*al*c2 
a(3,11  )*a (2 , 1 0 ) 

«(5,7)»a(2,10) 
a( 6 , 8  )«a ( 2 , 1 0 ) 

a(2,11 )»al*l/b4*(-»*c1e(b4+a“2)*e3) 

*(5,8 )*a( 2,11) 

a(3,2)-b4/l*c3 

a(6,5)-a(3.2) 

a(3 ,3)*c0-t*c2 

a(6,6)«a(3,3) 

a(3,10)«al/l*(c1-t*c3) 

«(6,7  )«a ( 3 , 10 ) 
a ( 7 , 5 ) ■ b4*  c2/al 

a(8 ,6  )«a(7,5 ) 

a(10,2)-a(7,5) 

a(11,3)*a(7,5) 

a(7,6)xsl/al*(-t*c1  +  (b4-ft“2)#c3) 

a ( 10 ,3  )**( 7 ,6 ) 

a(7,7)*c0-t*c2 

a(10,10)«a(7,7) 

a(7,8)xl*(c1-(s+t)*c3) 

a ( 1 0 , 1 1  )*a(7 ,8  ) 

a(8,5)-b4*(e1-a*c3)/al/l 

•  (11  ,2)»a(8,5) 
a(8,7)«b4*c3/l 
a ( 1 1 , 10  )*a ( 8 ,7 ) 
a( 8,8)«c0-a*c2 

•  (11 ,11  )*a(6 ,6  ) 

a(9,4)*-«u*l*(w*2)*sir(th)/th 

a ( 12 , 1  )*a( 9 .4 ) 

return 
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Table  K3.  Computed  Values  d  the  Determinant  d  Transfer  Matrices  far  Uniform  Rod  and 
Timoshenko  Beam  for  Various  Frequencies  of  Vibration* 


Frequency 
(rad/  sec) 

Computed  Values  d 
Transfer 

the  Determinant  d 

Matrix  for 

Longitudinal  Vibration 
d  Uniform  Rod 

Flexural  Vibration 
d  Uniform  Timoshenko 
Beam 

10 

. 

1 

100 

1 

1 

1000 

1 

0.99999 

10000 

1 

•0.0625 

100,000 

1 

1.80144  xlO16 

1000,000 

_ 

1 

— 

«• 

ft 


■  *  Far  the  raj,  elastic  modulus  is  7.46 XlO10  P3  (103  XlO4  pi),  cross -sectioul  area  is  6X14x10  3  ml 

(9375x15*  in*;. mass  density  is  72  kg/m 3  (26X10 "4  Ibm  /in*),  ard  iengdi  is  2.43 Xl5!  m  (9.75  in). 

Fcr  the  Timcehenko  beam,  elastic  modulus  is  7.46xl010  Pa  (10.8xl01D  psi),  shear  modulus  is  2.75X1010  Pa 
1 4 £) Xl0“  ail  cross-sectional  area  is  6X)4xl0"5nt:  t9375xl0"2  in*),  second  moment  d  area  inertia  is 
435x10*0  m4  (1  098  XlO"5  in4),  mass  density  is  72  kg/m*  123x10"*  Ibm/in *),  rodiia  d  gsranon  is 
2.71  xlO"3  m  ,10838x10"' in ), and  length  is  2,43x10"  nt  (9.75  in) 

**  Value  exceeded  the  capacity  d  the  computer 
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APPENDIX  L 

SON- DIMENSION AUZED  FORMS  FOR  TRANSFER  MATRICES  OF  A 
3-BAY  PLANAR  LATTICE  STRUCTURE 


The  non-dimensianalized  forms  for  the  two  transfer  matrices  Xx  and  X;  far  the  analysis 
of  wave  propagation  and  vibration  in  a  3-bay  planar  structure  are  derived  in  this  appendix. 

Fran  eqn,  (39)  and  eqn.  (40),  the  two  transfer  matrices  X;  and  X;,  corresponding  to  a 
state  vector  Z  such  that 


Z  = 


d, 

dtt 

Pn 

Pi 


where 


d  = 


u 

~w 


crd  where  p  — 


'  1 
M 

v  r 

N 


and  where  t he  subscripts  denote  the  member  numbers,  are  given  by 


and  X2  = 


Ci  0  0  c2 
0  Cl  Cj  o 
0  C)  C4  0 
C3  0  0  Cy 

I  0 

0  I 

G<CilGx  -G£{lC  jG3 

G£<C;xGl  G2(C)-C<C{lCi)Gi 


I  is  the  identity  matrix,  0  is  the  null  matrix,  C2f  C2,  C3  and  C«  are  3x3  submatrices  and  Gj,  G2,  G 
and  G«  are  transformation  matrices  such  that 


cat  0  0  0 

0  c0-<r2  *[ci-(<r+r)c3] 

0  $?_  c0^rc2 


(  sin0 
EA  0 


[-  aci  +  (p+o^cj]  0 
&EJ  0 


17  (Cl_TC3)  17 Cl 


~Ci 


-J-  t-r  ci  + 


■y^ 


g4^ 

e 

sinQ  o 


r  i  3*£J 

[Ci-CTCj]  C: 


Ci  =  i 


*[c!-(<m)c3]  o 

Cn-OT,  0 


Cj  Ca~<rC2 


VV^+I-m2  4(trh) 

,\(\|  cosh  Xi  +  XiJ  cos  Xj) 

j j 

A(—  sinh  \  sin  Xz) 


A(ccsh  Xi  -  cos  Xj) 
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^  GM 


_  pf2!2^ 
T  “  £J 


3*  =  \il— 

P  P  EJ 


Now  assume  that  the  connecting  bars  have  identical  square  crass -sections  cf  side 


that 


» 


S 


r 


-7 

04) 

.  _  a 

1  ~  M2 

(L5) 

A-?a2 

(L6) 

Fran  the  elastic  theory  of  isotropic  materials, 


G  = 


£ 

2(1-*) 


(L7) 


where  v  is  Poisson’s  ratio.  Now  let 


n!  ^aiSsL  (u) 

and 

€a  «  —  (L9) 

a 

where  fl  and  £a  are  both  non-diraensianalized  lumped  parameters.  Using  Equation  (L4)  through 
Equation  (L9),  the  variables  0,  c r,  t,  and  (3*  used  in  the  transfer  practices  Xx  and  X:  can  be 


rewritten  as 


and  where 


Equations  (L2)  and  (L3)  can  now  be  rewritten  as 
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where 


= 


and  X2  = 


C:  0_  0_  C: 
0  C,  C2  0 
0_  C3  Ct  o_ 
Ci  o  0  c< 


1  0 

0  1  _ 

G4C2G1  “ G4C2^CiGj 

G&C^Gi  G:(T3-CtC?C})G3 


0  0 
0  0 
;  0 
0  / 


Ct  = 


cos  Cl  0  0 

0  co  -or 2  [cj—  (£+t)cj] 

0  p*cj  c.  -  t c2 


Cz  = 


0 

Ci 

Cl  -  tc3 


■[aeJ+(0‘+oz)c3] 


sinQ 

n 

0 

0 


Ci 


C3  = 


C4  = 


0  [-TCi  +  (0‘+T2)c3] 

0  ^(Ct-OCj)  g*C2 

fisinfl  0  0 


Cq-^Ci  [Ci-(o»-t)C3]  0 

$C3  Cg  -  <r:  0 

0  0  CCS  fl 


(US) 


(L16) 


Notice  that  Xx  and  X2  in  Equations  (L15)  and  (L16)  correspond  to  a  non- 
dime  nsi  ana  lized  state  vector  Z  and  the  elements  in  the  two  matrices  are  functions  of  non- 
dimensionalized  lumped  parameters  l a  and  fl . 


N  .\ 
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